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Preface 


This  technical  report  is  concerned  with  the  theory 
and  observations  of  acoustic  wave  propagation  through  the 
atmospheric  boundary  layer.  The  observations  consist  of 
amplitude  and  doppler  shift  variations  recorded  at  a 
single  receiver  about  10  kilometers  from  the  transmitting 
site.  These  observations  provide  indirect  evidence  of 
the  magnitude  and  scale  of  the  wind  velocity  fluctuations, 
which  are  the  main  causes  of  the  acoustic  propagation 
variations.  Better  resolution  of  the  scale  of  these 
fluctuations  and  their  distributions  as  a  function  of 
height  can  be  achieved  by  incorporating  an  array  of 
receivers  recording  simultaneously,  and  we  are  presently 
working  at  this  extension.  The  array  data  can  also 
define  the  diurnal  variations  of  the  wind  and  temperature 
profiles. 
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ABSTRACT 

An  infra-sonic  wave  probe  has  been  constructed  to  in¬ 
vestigate  the  variability  of  the  wind  and  temperature  struc¬ 
tures  in  the  lower  troposphere.  The  probe  could  be  used  for 
detecting  Internal  atmopsheric  gravity  waves  and  for  studying 
air  pollution  meteorology.  The  turbulence  and  diurnal  varia¬ 
tions  of  the  atmospheric  boundary  layer  were  detected  by  a 
prolonged  operation  of  the  probe  with  a  fixed  receiver  locat¬ 
ed  9200  m  to  the  SE  of  the  source.  At  periods  of  .5  to  8 
minutes,  the  wind  fluctuations  and  eddy  sizes  are  inferred 
from  doppler  shifts  and  amplitude  variations.  The  root-mean- 
square  horizontal  wind  fluctuations  are  .13  +  .06  m/s  at 
heights  of  about  192  m.  The  horizontal  scales  of  the  eddies 
range  from  200  m  to  1400  m,  while  the  vertical  scales  are 
99  ♦  28  m.  Therefore  the  eddies  appear  to  be  horizontally 
elongated  in  the  atmospheric  boundary  shears.  Signal  ampli¬ 
tude  variations  at  periods  of  .5  to  6  hours  as  well  as  diurnal 
variations  of  air  temperature  and  winds  suggest  6  fundamental 
effective  wind  profiles  in  the  atmospheric  boundary  layer. 

The  signal  source  of  the  probe  is  a  tube-resonator  operated 
at  a  constant  frequency  of  around  13„5  cps,  and  the  receiver 
is  a  Globe  microphone  with  a  phase  lock  amplifier. 
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GLOSSARY  OF  TECHNICAL  TERMS 

Some  technical  terms,  which  are  often  used  in  this  work 
and  may  not  he  familiar  to  readers,  are  listed  and  briefly 
explained.  The  number  following  each  term  denotes  the  page 
where  the  term  is  defined.  Definitions  for  commonly-used 
meteorological  terms  can  be  found  in  Glossary  of  Meteorology 
(Huschke  1959). 

EFFECTIVE  WIND  is  defined  by 

ve  •  (vy  -  vyo)  +  (c  -  ce), 

where  Vy  =  the  horizontal  wind  component  along  the  source- 
receiver  line, 

0  =  the  speed  of  sound, 

and  the  subscript  "o"  denotes  the  value  at  the  Earth's 
surface.  (41) 

EFFECTIVE  WIND  SHEAR  is  the  vertical  gradient  of  the 
effective  wind. 

elevated  EFFECTIVE  WIND  SHEAR  is  the  positive  effective 
wind  shear  of  the  upper  layer  in  a  two-layer  model.  The  ef¬ 
fective  wind  shear  in  the  lower  layer  may  be  positive,  ne¬ 
gative,  or  zero,  (45) 

AVERAGE  RAY  HEIGHT  (Hav)  is  the  average  height  to  which 
the  received  signal  has  propagated.  (41) 

SHEAR  VANISHING  HEIGHT  (zm)  is  the  height  where  the 
positive  effective  wind  shear  vanishes.  The  shear  vanishing 


height  is  at  about  200  m  to  600  m  above  the  Earth’s  surface. 
(43) 


SIGNAL  PENETRATING  HEIGHT  (H)  is  the  maximum  height 
to  which  the  received  signal  has  propagated.  (41) 

FUidSING  FACTOR  (f)  is  the  ratio  of  the  observed 
signal  amplitude  to  the  amplitude  which  one  would  expect  if 
the  signal  was  propagating  over  the  same  horizontal  distance 
in  a  homogeneous  atmosphere.  The  focusing  factor  is  a 

measure  of  the  geometrical  spreading  of  neighboring  rays. 
(31)  (42) 
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LIST  OF  SYMBOLS 

(The  number  following  each  symbol  denotes  the  page 
where  the  symbol  is  defined  or  first  occurs.) 
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CHAPTER  1.  INTRODUCTION 

This  work  concerns  the  development  of  an  infra-sonic 
wave  probe  for  investigating  meso-scale  wind  and  temperature 
variations  in  the  lower  troposphere.  Wind  and  temperature 
fluctuations  have  been  measured  by  research  aircraft  and 
meteorological  towers,  but  these  direct  measurements  are 
limited  in  either  time  or  space.  For  continuously  monitoring 
average  turbulence  properties  over  a  large  area,  a  wave-pro¬ 
pagation  probe  seems  to  be  more  suitable.  When  strong  near¬ 
surface  wind  shears  or  inversion  layers  are  present,  one  may 
send  an  acoustic  signal  up  and  receive  it  some  distance  away 
in  the  downwind  direction.  The  amplitude  variations  and 
doppler  shifts  of  the  received  signal  will  give  information 
about  turbulence  in  the  atmospheric  boundary  layer,  when 
strong  jet  streams  exist  in  the  upper  troposphere,  one  might 
also  get  longer  range  signal  transmissions  by  using  the  large 
wind  shear  below  the  jet  stream  cores.  Then  the  signal  varia¬ 
tions  would  give  information  about  perturbing  winds  due  to 
internal  atmospheric  gravity  waves,  which  have  often  been 
found  to  accompany  atmospheric  jet  streams. 

A  quarter  wave  tube-resonator  with  resonant  frequency 
of  around  13.5  cps  was  built  as  an  acoustic  source.  A  Globe 
microphone  and  a  phase  lock  amplifier  were  used  as  the  receiver. 
Field  data  was  collected  with  signal  travel  distance  only  up 
to  9.2  km.  The  signal  penetrating  heiaht  was  limited  by  a 
shear  vanishing  height  which  is  about  200  m  to  600  m  above 


the  Earth's  surface.  From  the  observed  average  signal  ampli¬ 
tude  and  weather  information,  the  average  ray  height  was  es¬ 
timated  to  be  about  192  m. 

Signal  variations  with  periods  of  about  .5  to  8  min. 
were  mainly  caused  by  the  turbulent  air  motion  drifting  with 
the  average  wind  in  the  atmospheric  boundary  layer.  Wind 
perturbations  due  to  gravity  waves  were  ruled  out,  because 
the  few  gravity  waves  which  occurred  during  the  measurements 
did  not  seem  to  have  large  enough  fluctuating  wind  components 
along  the  source-receiver  line  of  the  probe  to  cause  detectable 
signal  variations.  A  ray  theory  for  signal  propagation  in  the 
atmospheric  boundary  layer  was  developed  to  analyze  the  field 
data.  Formulas  for  root-mean-square  amplitude  anc  phase  varia¬ 
tions  of  the  probing  signal  were  derived.  Then,  from  the  ob¬ 
served  periodicity  and  the  average  wind  component  along  the 
source-receiver  line,  one  estimated  the  horizontal  turbulence 
scale.  From  the  observed  phase  shift  and  the  estimated  hori¬ 
zontal  turbulence  scale,  one  inferred  the  average  horizontal 
wind  fluctuation.  Finally,  from  both  the  phase  shift  and  the 
amplitude  variation,  one  inferred  the  vertical  turbulence 
scale.  The  inferred  wind  fluctuations  and  turbulence  scales 
are  generally  consistent  with  those  obtained  from  meteorolo¬ 
gical  tower  measurements  elsewhere.  The  horizontally  elon¬ 
gated  shape  of  the  inferred  eddies  appears  to  be  characteris¬ 
tic  of  the  turbulence  in  boundary  shears. 

Signal  amplitude  variations  with  periods  of  .5  to  6 
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hours  suggest  diurnal  properties  of  the  atmospheric  boundary 
layer,  whose  wind  and  temperature  structures  can  be  described 
by  6  fundamental  effective  wind  profiles. 

Chapter  2  describes  the  equipment.  This  includes,  in 
addition  to  the  signal  source  and  receiver,  wind  screens  and 
a  Daniels  pipe  for  improving  the  sianal-to-noise  ratio  at  the 
receiver.  Chapter  3  presents  the  theoretical  aspect  of  the 
probe  experiment.  It  starts  with  a  summary  of  the  general 
ray  theory.  Then  it  introduces  the  shallow  angle  approxi¬ 
mation,  which  is  proper  for  signal  propagation  in  the  atmos¬ 
pheric  boundary  layer.  This  approximation  helps  considerably 
in  evaluating  the  ray  characteristics  of  6  fundamental  effec¬ 
tive  wind  profiles.  A  linear  theory  of  parabolic  profile 
fluctuations  is  also  presented  to  predict  possible  amplitude 
variations  and  doppler  shifts.  Finally,  in  chapter  3,  root- 
mean-square  amplitude  and.  phase  variations  due  to  boundary 
layer  turbulence  are  derived,  assuming  a  three-dimensional 
Gaussian  correlation  function.  Chapter  4  gives  experimental 
results.  Firstly,  typical  examples  of  field  data  and  its 
analysis  are  demonstrated.  The  accuracy  of  all  analyzed  data 
is  estimated.  Then  wind  fluctuations  and  turbulence  scales 
inferred  from  all  available  data  are  discussed  in  the  light 
of  current  status  of  the  research  about  atmospheric  turbulence 
The  possible  signal  variations  due  to  gravity  waves  are  dis¬ 
cussed.  Experimental  evidences  of  diurnal  variations  of  the 
atmospheric  boundary  layer  are  also  presented.  Chapter  5 
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summarizes  the  resuits  of  the  prote  experiment  and  suggests 

directions  for  further  Improvements.  The  appendices  include 

the  details  of  experimental  results  and  some  algebras  of  the 
theory. 


JhLiFIER  2.  THE  lNFR,v-3CNI2  W^VE  FR03E 


S 


In  order  to  use  the  amplitude  and  frequency  information 
of  the  probe,  one  must  have  a  signal  source  whose  power  output 
and  frequency  are  steady.  The  signal  receiver  must  be  sensi¬ 
tive  enough  to  measure  the  amplitude  and  frequency  variations 
due  to  the  atmospheric  fluctuations  that  are  to  be  investi¬ 
gated.  The  signal  source  was  obtained  by  building  a  quarter 
wave  tube-resonator  and  closely  controlling  the  driving  fre¬ 
quency.  The  signal  receiver  consisted  of  a  Globe  microphone 
and  a  phase  lock  amplifier  which  can  measure  variations  of 
the  acoustic  amplitude  down  to  .001  jih  (dvne/cm^)  and  doppler 
shifts  accurate  to  .01  i  of  the  source  frequency.  The  perfor¬ 
mance  of  the  signal  receiver  very  much  depends  on  the  signal- 
to-noise  ratio  at  the  receiving  site.  With  a  fixed  source 
power,  an  improvement  of  signal-to-noise  ratio  can  only  be 
obtained  by  reducing  the  background  noise,  wind  screens  and 
a  Laniels  pipe  were  found  very  helpful  in  this  respect. 
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2-1.  '."he  signal  source 

The  signal  is  produced  by  a  quarter  wavelength  reso¬ 
nant  tube,  which  is  a  Sontube  20  feet  long  and  20  inches  in 
diameter,  sittino-  upright  on  a  base  frame  3  feet,  above  the 
ground.  The  upper  end  of  the  tube  is  open  and  its  bottom 
end  is  closed  by  a  rubber  diaphragm.  Two  aluminum  plates, 
one  on  each  side  of  the  diaphragm,  are  pinched  together  as 
a  piston.  The  piston  is  driven  up  and  down  by  a  one  horse¬ 
power  L-2  motor  through  a  linkage,  which  transfers  the  rota¬ 
tion  of  the  motor  shaft  into  the  vertical  motion  of  the  pis¬ 
ton. 

The  apparent  source  power  as  determined  by  a  receiver 
on  the  ground  is  4  times  the  real  source  power  because  of 
the  ground  reflection  which  doubles  the  sianal  amplitude  at 
the  receiver.  If  the  quarter  wavelength  resonant  tube  was 
situated  on  a  flat  plain  in  a  uniform  atmosphere  ,  the  ground 
reflection  would  make  the  signal  source  look  like  a  cipole 
source  with  a  separation  of  a  half  wavelength  between  the 
real  point  source  and  its  image  (horse  and  Ingard  196f ,  p. 
36c).  The  apparent  source  power  to  a  receiver  above  the 
crounc  would  be  4  T  cos^  [(7T/2)  sin  o(0]  >  where  ?  is  the 
real  source  power  emitted  from  the  top  of  the  tube,  and  o(c 
is  the  elevation  angle  of  the  receiver  in  a  coordinate  sys¬ 
tem  with  the  tube  bottom  as  the  origin.  Therefore,  to  a 
receiver  on  the  --round  { o(Q  =  0),  the  apparent  source  cower 
is  4  P.  Throuohout.  this  exoeriment ,  the  receiver  is  always 
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located  on  the  ground. 

At  resonance,  the  real  acoustic  power  fror  -pan 

enc?  of  the  tube  was  determined  to  be  15  watts.  1  he  quality 
factor  (3)  of  the  resonator,  defined  as  the  ratio  of  the 
resonant  frequency  to  the  frequency  width  of  half -power 
points  on  the  power  versus  frequency  curve  (Fig.  2-1.1),  was 
found  to  be  30.  _ne  can  analyze  the  resonant  tube  by  assuming 
two  plane  waves  inside  the  tube,  one  going  upward  and  the 
other  downward,  '.he  lower  boundary  condition  is  that  the 
particle  velocity  is  equal  to  the  piston  velocity,  the  upper 
boundary  condition  is  that  the  ratio  of  acoustic  pressure  to 
particle  velocity  is  equal  to  the  radiation  impedance  as 
given  by  Morse  (1946,  p.  333)  or  Morse  *  Ingard  (1966,  p. 

473).  Ihe  results  give 

Q  (the  quality  factor) 

=  (lT/4)pc/(.ia  +  ,qb), 

P0  (the  open  end  acoustic  power  at  resonance) 

=  (TT/2)  ?.a  ap  up)2  /  (<a  +  ib)2 

=  (8/IT)  Ra  (ap  up  3)2, 

and  P  (the  piston  end  acoustic  power  at  resonance) 
w 

=  nr/2)  (pc  ap  up)2  /  na  + 

=  2  P  C  3  (ap  up)2, 

where  p.  =  the  characteristic  impedance 
=  407  newton-sec/m  , 

?.a  =  the  radiation  resistance  at  the  open  end 


Fig.  2- 


.1  The  relative  acoustic  power  vs. 

frequency  of  the  signal  transmitter 


SOU  Relative 
I  acoustic 
I  power 
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=  d/2)  e  z  (Ka)2 

3 

=  .829  newton-sec/m  , 
k  =  2  TT/\, 

=  the  acoustic  wavelength 
=  25  m, 

a  =  the  tube  radius 
=  .254  m, 

ap  =  the  piston  radius 
=  . 238  m. 

Up  =  the  piston  velocity 
-  .269  m/s, 

?.b  =  r  l  (0  c, 

=  the  dissipation  resistance  inside  the  tube  due 
to  air  leakaqes,  nonlinear  effects,  eddy  and 
molecular  viscosity,  thermal  conduction,  wall 
flexibility,  etc. 

7 f  =  the  imaginary  part  of  the  propagation  constant 
inside  the  tube, 

L  =  the  effective  tube  length 

=  A/4. 

The  calculation  results  can  be  illustrated  by  the  numbers 
listed  in  Table  2-1.1.  The  above  theory  does  not  explain 
the  observed  3  and  power  output  of  the  transmitter.  For  an 
observed  acoustic  power  output  of  15  watts,  the  Q  would  be 
42  instead  of  the  observed  30,  which  corresponds  to  an  out¬ 
put  power  of  only  7.F  watts.  The  possible  reasons  for  the 
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Table  2-1.1  r-.he  predicted  quality  factor,  acoustic  power, 

and  acoustic  pressure  ot  the  sicnal  transmitter 


2 

P_,  (watt) 

O 

P  (watt) 

P  (bar) 
c 

25 

5.4 

84 

.035 

30 

• 

no 

100 

.042 

35 

ii 

120 

.050 

42 

15 

140 

.056 

386 

1290 

1290 

.54 

Remarks : 

2  =  the  quality  factor 

PQ  =  the  open  end  acoustic  power  in  resonance 
Pp  =  the  piston  end  acoustic  power  in  resonance 
Pp  =  the  piston  end  acoustic  pressure  in  resonance 

discrepancy  between  the  theory  and  observations  are; 

(1)  The  actual  radiation  resistance  could  be  greater 
than  what  is  assumed  in  the  theory,  which  is  derived  for 
a  flanged  piston. 

(2)  diome  unaccounted  process  in  the  system  could  have 
lowered  the  2,  without  simultaneously  lowerina  the  power 
output. 

(3)  The  sensitivity  of  the  receiver  microphone  could 
be  inaccurate  so  that  the  signal  amplitude  has  been 
overestimated . 

If  there  was  no  dissipation  inside  the  tube,  the  .t  would  he 

386,  the  resonant  oower  output  1290  watts,  anc  the  acoustic 


pressure  .54  bar,  according  to  the  theory.,  'ihis  seems  to  be 
not  only  impossible  because  of  the  large  acoustic  pressure 
(.54  bar)  which  could  have  caused  nonlinear  dissipation,  but 
also  impractical  because  of  the  laroe  j.  .^ith  a  j  of  3 Pf, 
the  output  power  would  decrease  by  504  when  the  transmitter 
frequency  deviates  from  the  resonant  frequency  by  .  13»,  which 
could  easily  be  effected  by  an  air  temperature  variation  of 
only  .  P°J. 

Iwo  control  modes  of  the  signal  transmitter  are  now 
available,  i.e.,  resonance  control  mode,  and  frequency  con¬ 
trol  node.  In  the  resonance  control  mode,  the  motor  speed 
is  fixed  at  resonance  by  controlling  the  piston  velocity  to 
be  always  in  phase  with  the  acoustic  pressure  at  the  bottom 
of  the  tube.  Since  the  resonant  frequency  varies  with  the 
air  temperature,  the  signal  source  under  the  resonance  control 
mode  does  not  necessarily  have  a  constant  frequency,  and, 
as  a  result,  is  not  suitable  for  studying  doppler  shifts  pro¬ 
duced  by  propagation  in  the  atmosphere.  In  the  frequency 
control  mode,  the  motor  speed  is  controlled  to  run  at  a  con¬ 
stant  frequency,  which  can  be  readjusted  to  be  equal  to  the 
resonant  frequency  as  the  air  temperature  changes,  tfith  a 
.j  of  30,  the  source  power  can  often  be  kept  within  10,-i  of 
the  resonant  power  by  readjusting  the  motor  speed  once  to 
thrice  over  a  whole  night.  Even  during  periods  with  maxi¬ 
mum  time  rates  of  change  of  the  air  temperature  on  a  clear 
summer  day,  i.e.,  about  +  3°J  per  hour,  the  source  power  can 
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still  be  kept  within  15 o  of  the  resonant  power  by  readjusting 
the  motor  speed  once  per  hour.  In  this  experiment,  the  fre¬ 
quency  control  mode  was  used  almost  exclusively.  The  elec¬ 
tronic  systems  for  the  motor  control  at  the  transmitter  and 
for  the  data  processing  at  the  receiver  are  very  similar, 
both  using  phase  lock  amplifiers.  A  schematic  diagram  (Fig. 
2-2.4)  in  the  next  section  shows  the  major  constituents  of 
the  entire  probe.  Letails  of  the  transmitter  control  should 
always  be  tried  out  to  fit  the  special  purpose,  for  which 
the  received  signal  is  processed.  For  example,  in  the  last 
stage  of  this  experiment,  the  received  doppler  shifts  were 
integrated  to  directly  record  signal  phase  shifts  originated 
by  atmospheric  turbulence  (dection  2-2).  To  cesicn  the  pro¬ 
per  transmitter  control  for  such  an  operation  should  minimize 
deviations  of  the  integrated  frequency  error. 

The  transmitter  frequency  can  usually  be  controlled  to 
stay  within  .005 i  of  the  average  frequency  (Figs  2-1.2  anc 
2-1.3).  Fig.  2-1.2  compares  the  source  frequency  and  ampli¬ 
tude  with  the  sicnal  frequency  and  amplitude  received  on  the 
around,  about  20  m  away  from  the  transmitter.  The  source 
frequency  and  amplitude  were  monitored  by  suspending  a  ceramic 
microphone  near  the  piston  end  of  the  tube  resonator.  The 
sensitivity  of  the  ceramic  microphone  has  not  yet  been  cali¬ 
brated.  Therefore,  the  source  amplitude  is  expressed  in  an 
arbitrary  scale.  In  rig.  2-1.2,  the  frequency  records  of  the 
transmitter  and  the  receiver  look  almost  identical,  with  the 
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receiver  so  close  to  the  source.  Large  deviations  of  the 
transmitter  frequency  can  be  caused  by  the  wind  blowing  over 
the  mouth  of  the  acoustic  resonator,  electric  surges  in  the 
power  line  to  the  transmitter  control  system,  etc.  Therefore, 
to  interpret  the  received  doppler  (or  phase)  shifts,  one 
should  always  first  compare  them  with  the  source  frequency 
record  and  count  out  any  variations  which  are  due  to  the 
source.  Fig.  2-1.3  compares  the  source  frequency  with  the 
signal  frequency  and  amplitude,  which  were  received  9200  m 
to  the  SE  of  the  source.  In  order  to  determine  the  back¬ 
ground  noise,  the  transmitter  was  scheduled  to  be  on  for  22 
min.  and  off  for  2  min.,  with  a  cycle  per  24  min.  After  being 
corrected  for  source  frequency  deviations,  the  received  fre¬ 
quency  variations  can  be  either  due  to  medium  fluctuations  or 
due  to  the  background  noise.  Reliable  doppler  shifts,  which 
are  caused  by  wind  and  temperature  fluctuations,  can  only  be 
determined  after  one  fully  understands  noise  effects  on  the 


receiver  behavior 
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2-2.  The  signal  receiver 

The  signal  receiver  includes  three  major  components s 
a  Globe  microphone,  a  filter-amplifier  and  a  phase  lock  sys¬ 
tem. 

The  Globe  microphone  is  a  sensitive  capacitor  micro¬ 
phone.  Its  frequency  response  is  flat  within  3  db  from  .1 
cps  to  450  cps.  It  has  a  sensitivity  of  .0225  v/ yjb.  The 
microphone  output  is  sent  to  the  filter-amplifier. 

The  filter-amplifier  had  a  4  cps-bandwidth  centered  at 
13.5  cps.  Later  the  bandwidth  was  narrowed  to  1  cps.  Its 
amplification  gain  ranges  from  10  to  5000.  The  filter  output 
is  rectified  and  averaged.  This  output  is  denoted  as  the  A. 2, 
level,  which  represents  the  total  level  of  the  signal  and  the 
noise  within  the  4  cps-  or  1  cps-bandwidth.  The  filter  out¬ 
put  is  also  sent  to  the  phase  lock  system. 

The  phase  lock  system  phase-detects  the  signal  against 
a  relaxation  oscillator.  The  phase  error  is  integrated  to 
produce  a  frequency  error  signal.  This  error  signal  is  used 
to  keep  the  relaxation  oscillator  in  step  with  the  signal  and 
is  also  used  as  a  monitor  of  the  signal  frequency.  The  re¬ 
laxation  oscillator  also  operates  a  synchronous  detector 
through  which  the  audio  signal  is  passed.  The  output  of  this 
detector  is  averaged  and  denoted  as  the  synchronous  output. 

In  a  later  stage,  the  error  signal  was  integrated  to  give 
the  phase  shift  of  the  received  signal.  The  theoretical 
characteristics  of  the  receiver  phase  lock  loop  are  the 
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following  (Gardner  1967): 

The  effective  bandwidth  of  synchronous  output 
=  .05  cps. 

The  amplitude  of  synchronous  output 
=  (signal  amplitude)  ^cos  <e  noised  ' 
where  (0noise)  *  the  signal  phase  jitter  due  to  noise. 

The  mean  square  phase  jitter  due  to  noise 

^ noise ^ 

“  ^  noise )  /  ^0  psignal^' 
where  (Pnoise^  =  noise  power, 

(psignal^  =  the  si9nal  power. 

The  maximum  doppler  shift  rate  that  can  be  tracked  when 

^signal )  >  *pnoise^ 

=  .012  cps/sec. 

The  minimum  signal-to-noise  ratio  for  lock 
=  1/5. 

The  receiver  has  been  constantly  under  improvement. 

For  the  experimental  results  presented  in  this  work,  one  can 
distingnish  three  stages  of  the  receiver  condition: 

1)  From  May  8  to  Nov.  16  in  1970.  the  receiver  had 
a  wide  band  (4  cps)  filter  and  recorded  frequency  shifts. 

2)  From  Nov.  16  to  Nov.  24  in  1970,  the  receiver 
had  a  narrow  band  (1  cps)  filter  and  still  recorded 
frequency  shifts. 

3)  From  Nov.  24,  1970  to  Jan.  7,  1971,  the  re¬ 
ceiver  had  a  narrow  band  filter  and  recorded  phase  shifts. 
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The  receiver  characteristics  are  demonstrated  in  Figs. 
2-2.1  to  2-2.3  for  the  three  different  stages.  For  signal- 
to-noise  ratios  of  about  1,  the  frequency  is  steady  over 
periods  of  more  than  1  min.,  although  over  periods  of  less 
than  1  min.,  there  are  frequency  fluctuations  due  to  the  noise. 
Therefore,  the  minimum  signal-to-noise  ratio  for  detecting 
frequency  (or  phase)  shifts  with  periods  longer  than  1  min. 
is  about  unity  for  all  three  stages. 

The  entire  probe  system  is  schematically  shown  in  Fig. 

2-2.4. 
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Fig.  2-2.1  Receiver  characteristics 

(4  cps-bandwidth  and  doppler  shifts) 
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Fig.  2-2.2  Receiver  characteristics 

(1  cps-bandwidth  and  doppler  shifts) 
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Fig.  2-2.3  Receiver  characteristics 

(1  cps-bandwidth  and  phase  shifts) 
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Fig.  2-2.4  Flow  diagram  of  the  probe 
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2-3.  Wind  screens  and  the  Daniels  pipe  for  improving  the 

signal- to-noise  ratio  at  the  receiver 

The  improvement  of  the  signal-to-noise  ratio  at  the 
receiver  is  imperative  in  carrying  out  the  experiment  success¬ 
fully.  In  addition  to  the  electronic  phase  lock  technique, 
which  proves  to  be  unique  in  tracking  the  signal,  some  me¬ 
chanical  devices  are  still  necessary  to  reduce  the  background 
noise  at  the  receiver.  There  are  two  kinds  of  background 
noises  the  random  noise  and  the  acoustic  noise.  The  random 
noise  is  mainly  caused  by  the  local  wind.  The  acoustic  noise 
is  mainly  caused  by  the  high  speed  automobile  traffic.  Wind 
sreens  made  of  silk  can  reduce  the  wind  noise  by  inhibiting 
the  air  flow.  A  Daniels  pipe  (Daniels  1959)  with  wind  screens 
can  reduce  both  random  and  acoustic  noises. 

A  portable  wind  screen  made  of  two  cylindrical  silk 
layers,  one  outside  the  othfer,  has  been  used  for  exploratory 
field  trips.  The  outer  cylinder  is  2  feet  high  and  2.5  feet 
in  diameter,  and  the  inner  one  is  1.5  feet  high  and  about  1 
foot  in  diameter.  The  exact  gain  of  the  wind  screen  is  yet 
to  be  measured.  The  noise  reducing  capability  seems  to  depend 
on  the  wind  speed.  The  signal  is  unattenuated,  by  the  screen. 

The  Daniels  pipe,  which  acts  as  a  line  microphone,  is 
a  plastic  pipe  25  meters  long  with  an  inside  diameter  of  1 
inch.  One  end  of  the  pipe  is  fitted  to  the  Globe  microphone 
of  the  receiver,  and  the  other  end  pointed  at  the  signal 
source.  Evenly  distributed  along  the  pipe  are  25  leaks,  which 
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are  made  of  capillary  tubes  of  different  sizes*  The  leaks  at 
the  source  end  of  the  pipe  have  larger  diameters  and  smaller 
lengths  than  the  leaks  at  the  microphone  end,  so  that  the 
waves  coming  through  all  the  leaks  will  arrive  at  the  micro¬ 
phone  with  approximately  the  same  amplitude.  The  sound  speed 
in  the  pipe  is  equal  to  that  in  the  air.  Therefore,  for  the 
signal  travelling  in  the  direction  of  the  pipe,  each  leak 
admits  a  wave  which  is  in  phase  with  the  waves  from  all  the 
other  leaks.  The  wind  noise,  however,  because  of  its  limited 
spacial  coherence  is  incoherently  summed  in  the  pipe.  Also 
the  acoustic  noise  coming  from  other  directions  is  attenuated 
because  the  waves  picked  up  at  different  leaks  are  out  of 
phase.  The  theoretical  amplitude  response  as  a  function  of 
the  direction  is  shown  in  Fig.  2-3.1  (Clson  1947,  p.  280). 

The  Daniels  pipe  as  well  as  the  Globe  microphone  are  buried 
and  the  leaks  are  covered  with  small  wind  screens.  This  sys¬ 
tem  reduces  the  background  noise  by  a  factor  of  about  3.  The 
signal  is  essentially  unattenuated  by  the  Daniels  pipe. 


i  v  w  v 
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CHAPTER  3.  A  RAY  THEORY  FOR  SIGNAL  PROPAGATION  IN  THE 
ATMOSPHERIC  BOUNDARY  LAYER 

A  simple  ray  theory  is  developed  to  analyze  the  field 
data  of  the  experiment.  With  horizontal  travel  distances  of 
5  to  10  km,  the  signal  penetrating  height  is  limited  by  the 
atmospheric  boundary  shear,  which  is  200  to  600  m  thick. 

Because  of  the  small  elevation  angles  of  the  propagating  ray 
(0(^.3  rad.),  the  theory  adopts  a  shallow  angle  approximation. 
One-dimensional  effective  wind  profiles  are  proposed  to  inter¬ 
pret  long-period  signal  amplitude  variations  with  periods  of 
.5  to  6  hours.  From  these  models,  the  average  ray  trajectories 
and  penetrating  heights  can  be  estimated.  Superposed  on  the 
long-period  atmospheric  variations  are  short-period  wind  and 
temperature  fluctuations  with  periods  of  .5  to  8  min.,  which 
cause  short-period  amplitude  variations  and  phase  (or  doppler) 
shifts.  Horizontal  wind  fluctuations  and  turbulence  scales 
can  be  inferred  from  the  short-period  signal  variations  by  a 
three-dimensional  fluctuation  model. 

The  first  section  of  this  chapter  (Section  3-1)  sum¬ 
marizes  the  general  ray  theory,  whose  details  can  be  found  in 
Hayes  (1970)  and  Lighthill  (1965).  Section  3-1  also  discusses 
various  signal  attenuating  processes  in  the  atmosphere,  dif¬ 
ferent  ways  of  evaluating  ray  integrals,  and  finally  the  multi- 
path  ray  interference.  Section  3-2  introduces  one-dimensional 
atmospheric  models  and  the  shallow  angle  approximation.  These 
help  simplify  the  algebra  of  integrations.  Section  3-3 


presents  fundamental  effective  wind  profiles  which  are  sug- 
gested  from  the  long-period  experimental  results  (section  4-4). 
Section  3-3  also  summarizes  the  ray  characteristics  of  these 
fundamental  profiles.  Section  3-4  is  a  linear  theory  of  pro¬ 
file  fluctuations  which  illustrates  the  relation  between  am¬ 
plitude  variations  and  doppler  shifts.  Section  3-5  deals 
with  signal  variations  due  to  short-period  wind  fluctuations. 
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3-1.  Ray  theory 

The  important  results  of  general  ray  theory  can  be  sum¬ 
marized  in  the  following  three  statements. 

a)  The  group  velocity,  which  describes  the  ray  trajec¬ 
tory,  is  the  gradient  of  frequency  with  respect  to  wave  vector 
in  the  augmented  space  (Hayes  1970).  The  augmented  space, 
which  comprises  wave  vector  Jc,  position  vector  t,  and  time  t, 
is  distinguished  from  the  propagation  space  (r,  t).  The  fre¬ 
quency  in  the  augmented  space  is  denoted  byj?(Jc,  r,  t),  and 
that  in  the  propagation  space  by  6U  ( ?,  t).  The  dispersion 
relation  is  then 

t),  (la) 

which,  for  the  acoustic  wave  in  a  moving  medium,  is 

60  =  kC  +  k-V,  (lb) 

where  C  =  sound  speed, 

V  =  wind, 

and  both  C  and  ^  are  functions  of  r  and  t.  Therefore  the 
acoustic  group  velocity  is 

=  nC  +  V,  (2) 

where  ft  =  unit  wave  vector. 

b)  The  time  derivative  of  wave  vector  along  the  ray  is 
the  negative  gradient  of  frequency  with  respect  to  position 
vector  in  the  augmented  space,  i.e.. 


sat  =  _n 

dt  6r» 


(3) 


where  ^_  =  2—  +  —  .  T7 

At-  ‘JM-  dt  V . 


dt  Vt  '  dt  '  v  r* 
c)  The  time  derivative  of  frequency  along  the  ray  is 
the  partial  derivative  of  frequency  with  respect  to  time  in 
the  augmented  space ,  i.e.. 


dO)  -  !£  ♦  $  .  2$ 

dt  “  *  *«t  *  ^t  * 


(4) 


The  coordinate  system  used  for  the  probe  is  illustrated 
in  Pig.  3-1.1.  The  wave  vector  magnitude  and  components  are, 

k  -  _ QL _  (5) 

C  +  Vx  cos  d  cos  /3  +  Vy  COs  of  sin  p  +  Vz  sin  d ' 

kx  =  k  cos  ol  cos  (3  # 

*  =  k  COs  ol  sin  (3  , 

kz  =  k  sinol  . 

The  signal  amplitude  is  calculated  by  considering  the 
geometrical  spreading  of  the  rays.  The  doppler  shift,  which 
is  usually  less  than  . 1 %  of  the  source  frequency,  can  be  ne¬ 
glected  in  evaluating  ray  trajectories.  The  apparent  source 
power  is  4  times  the  real  power  output  ,p,  of  the  signal  trans¬ 
mitter  because  of  the  ground  reflection  at  the  receiver 
(Section  2-1).  Then  the  power  emitted  in  the  solid  angle 
element  cos  dQ  d  0lo  d  ^3  0  is  (P/7T)  cos  ctQ  d  0(o  d  (30#  which, 
according  to  the  law  of  energy  conservation,  should  be  equal 
to  [p2/(2f  C)]  |n  •  x  Jd  d(30  at  tne  receiver/ 

where  p  is  the  peak  amplitude  of  acoustic  pressure, P  the  air 

density,  and  I  n  •  x  |  <a of  dA_  the  cross-section  of 
1  Wo  "dpo  I  0,0 


O*  ft.  N  <  X  »  W 

*  m  c\ 
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the  ray  tube.  Henqe  the  signal  peak  amplitude  is 

P  =  (2peocoS£/o/(n|*-||-x^-|  >)\  <6) 

A  convenient  way  to  express  the  signal  amplitude  received 
away  from  the  source  is  to  specify  its  focusing  factor,  which 
is  defined  as  the  ratio  of  the  actual  amplitude  to  the  ampli¬ 
tude  one  would  receive  if  the  medium  was  uniform  and  the 
signal  could  propagate  along  a  straight  line,  i.e.. 


f  =  R  /  I  n  •  x  2^-.  I 

I  Wa  Ma  I 


(7) 


«o  ^|3o 

The  signal  can  be  attenuated  by  other  processes  than 
the  geometrical  spreading,  which  are,  in  the  order  of  impor¬ 
tance  to  this  e^eriment:  vibrational  relaxation  of  atmos¬ 
pheric  molecules,  turbulence  scattering,  relaxation  effects 
of  fogs,  and  classical  absorption  due  to  viscous  forces  and 
heat  conduction.  The  amplitude  attenuation  caused  by  these 
processes  is  estimated  according  to  theoretical  studies  and 
laboratory  tests  reported  in  the  literature. 

The  major  signal-attenuating  relaxation  process  in  the 
atmosphere  is  due  to  the  vibrational  mode  of  oxygen  molecules 
(Kneser  1965).  The  water  vapor  content  decisively  controls 
the  absorption  as  shown  in  Fig.  3-1.2.  For  example,  the 
maximum  amplitude  loss  of  3  356  over  a  9.2  km  path  occurs  ap- 
proximatly  at  -10°C  with  a  relative  humidity  of  100%,  at  0°C 


with  a  reL  -ive  humidity  of  50%,  or  at  10°c  with  a  relative 


humidity  of  20%. 
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The  amplitude  loss  due  to  turbulence  scattering  depends 
on  the  intensity  of  boundary  turbulence,  which  can  be  esti¬ 
mated  from  the  3  m-level  wind  fluctuation.  The  estimated 
signal  amplitude  losses  for  three  turbulent  conditions  are 
listed  in  Table  3-1.1.  For  a  travel  distance  of  9.2  km,  the 
signal  amplitude  can  lose  11%  in  a  severe  turbulence  (i.e., 
with  a  rms  wind  fluctuation  of  1.5  m/s  at  3  m). 

The  signal  is  attenuated  in  fogs  because  of  an  irrever¬ 
sible  energy  transfer  from  the  signal  w-*vr  to  water  droplets 
and  the  saturated  air,  similar  to  the  relaxation  attenuation 
by  molecular  vibrational  modes.  Accox.^ng  to  a  theoretical 
study  by  Cole  and  Dobbins  (1970),  the  energy  attenuation  of 
waves  with  frequencies  less  than  130  cps  is  about  8  x  10  4 
neper/m  when  the  visibility  in  the  fog  is  24  m.  Therefore 
the  amplitude  loss  over  9.2  km  in  a  dense  fog  will  be  almost 
total.  This  might  have  actually  happened  to  the  probing 
signal  in  the  night  of  Oct.  9-10,  1970,  when  the  signal  am¬ 
plitude  received  at  the  9.2  km  site  gradually  vanished  as 
the  visibility  in  the  fog  decreased  to  about  24  m.  However 
the  loss  of  signal  that  night  could  also  be  interpreted  as 
due  to  the  loss  of  wind  shear  which  was  recorded.  The  rest 
of  the  field  data  analyzed  in  this  work  was  collected  without 
significant  fogginess. 

The  amplitude  loss  due  to  viscous  forces  and  heat  con¬ 
duction  of  the  air  is  estimated  to  be  about  only  .003%  over 
a  9.2  km  path  (Piercy  1969),  and  is,  therefore,  unimportant. 


Table  3-1.1  signal  amplitude  loss  due  to  turbulence 
scattering 


(Tv 

(m/s) 

.5 

1. 

wwm 

e 

02591 

.0013 

.0053 

.012 

Signal 

ampli¬ 

tude 

loss 

(Y  = 

4.5  km) 

296 

556 

996 

(Y  = 

9.2  km) 

456 

1096 

1796 

Remarks: 


1.  ~  the  rms  wind  fluctuation  at  3  m- level. 

2.  £  =  the  turbulent  energy  dissipation  rate. 

=  u*2|2 
3z  $ 

where  u*  =  the  friction  velocity 
~  Oy/2.5, 

_ 

ana  ^  =  the  vertical  wind  shear 

=  .033  sec""*-. 

3.  Y  =  the  travel  distance. 

4.  The  signal  amplitude  loss  is  estimated  by 
assuming  a  homogeneous  isotropic  turbulence 
with  outer  scale  of  200  m,  a  temperature  of 
10°C,  and  turbulent  energy  dissipation  rates 
estimated  above.  The  formulas  can  be  found 
in  Batchelor  (1957)  and  Tatarski  (1961). 


35 


To  evaluate  the  ray  trajectory  and  the  ray  parameters, 
e.g.,  £0  ,  and  ,  there  are  three  different  ways 

of  integration.  They  are  time-integration,  z-integration, 
and  y-integration. 

Time-integration,  as  described  by  Wesson  (1970),  is 
suitable  in  numerical  calculation. 

Z-integration  is  preferred  for  obtaining  analytical 
expressions,  when  atmospheric  parameters  are  assumed  to  be 
functions  of  z  only.  Since  is  now  not  an  explicit  func¬ 
tion  of  x  and  y,  one  gets  from  Eq.  (3) 

^  =  kox  =  constant,  (8a) 

and  ky  =  koy  =  constant,  (8b) 

for  each  ray.  Also  as  a  matter  of  fact,  the  doppler  shift 
is  usually  less  than  .1%  of  the  source  frequency.  Therefore, 
the  frequency  can  be  practically  considered  as  constant  when 
ray  trajectories  are  evaluated,  i.e., 

60  =  00Q  =  constant.  (9) 

From  Eqs.  (8)  and  (9),  the  Snell's  law  in  a  moving  medium  is 
obtained  as 


cos  o( 


vy  = 


cos 


(10) 


The  vertical  component  of  the  wind,  V2,  is  negligible  in  de¬ 
termining  the  ray  trajectory,  z-integration  is  used  in 
Sections  3-2  through  3-4. 

Y-integration  is  preferred  for  calculating  the  root- 
mean-square  value  of  signal  variations  due  to  random  medium 
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fluctuations  and  is,  therefore,  used  in  Section  3-5. 

One  should  beware  of  the  possiblity  of  receiving 
multi-path  rays.  For  example,  with  a  two-layer  atmosphere 
model,  one  often  predicts  two  rays  with  different  initial 
elevation  angles  landing  at  the  same  receiving  site  (Fig. 
3-3.4).  The  multi-path  rays  could  interfere  with  each  other. 
However,  from  the  experimental  results  (Chapter  4),  one  finds 
that  the  multi-path  ray  interference  is  negligible  in  this 
experiment.  There  are  two  reasons  for  thiss 

(a)  Firstly,  observed  phase  shifts  (in  radians)  are 
mostly  much  greater  than  observed  fractional  amplitude 
variations  (Tables  A-1..1  to  A-1.9),  while,  for  the  multi- 
path  ray  interference,  one  would  predict  phase  shifts  to 
be  only  slightly  greater  than  or  equal  to  fractional  ampli¬ 
tude  variation.  For  example,  one  considers  the  interference 
of  two  rays  with  amplitudes  respectively  of  1  and  r,  where 
r  is  smaller  than  1.  Their  phase  difference  is  assumed 
to  vary  randomly  between  0  and  2  7T  .  Then,  according  to 
vector  summation,  one  gets  the  rms  phase  shift  as 
0$  =  (r  /  2h)  (1  +  r2  /  8), 
and  the  rms  fractional  amplitude  variation  as 
0£/p  =  (  r  /  2h)  (1  -  r2  /  32). 

Hence,  for  the  double  ray  interference,  phase  shifts  are 
predicted  to  be  only  slightly  greater  than  fractional 
amplitude  variations.  /7hen  the  number  of  interfering  rays 
becomes  infinitely  large,  the  signal  variations  could 
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possibly  be  predicted  by  the  diffraction  theory  for  signal 
propagation  in  a  random  medium  (Tatarski  1961,  p.  185), 
and  phase  shifts  likely  tend  to  be  egual  to  fractional 
amplitude  variations.  However,  experimental  results 
give  phase  shifts  to  be  mostly  much  greater  than  frac¬ 
tional  amplitude  variations  (Fig.  3-1.3).  In  average, 
phase  shifts  are  measured  to  be  twice  as  large  as  frac¬ 
tional  ampliutde  variations.  Therefore,  the  multi-path 
ray  interference  seems  to  be  negligible  in  this  experi¬ 
ment. 

(b)  Secondly,  because  of  the  topography  of  the  experi¬ 
ment  site  (Fig.  4.1),  rays  reflected  from  the  ground  sur¬ 
face  at  closer  distances  often  restart  with  greater  ele¬ 
vation  angles  than  their  trapping  angles,  and  are  unlikely 
to  bend  down  again  so  as  to  interfere  with  the  direct  ray 
at  the  receiver. 

Therefore,  to  the  first  approximation,  one  can  neglect  the 
interference  of  multi-path  rays,  and  consider  only  one  direct 
ray. 

Due  to  atmospheric  inhomogeneities,  the  receiver  can 
also  have  scattered  rays,  which  interfere  with  the  direct  ray 
and  cause  diffraction.  Cne  will  see  in  Section  4-1.  that 
the  field  data  of  this  experiment  can  still  be  analyzed  by 
the  ray  theory  and  that  the  effect  of  diffraction  phenomena 
on  these  results  can  be  estimated.  The  horizontal  wind  fluc¬ 
tuation  inferred  with  the  ray  theory  is  possibly  smaller  than 
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the  true  value  by  33 %.  The  accuracy  of  inferring  the  ver¬ 
tical  turbulence  scale  with  the  ray  theory  is  generally  good 
except  for  a  few  cases  where  the  phase  shift  is  almost  equal 
to  the  fractional  amplitude  variation. 
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3-2.  One-dimensional  atmospheric  model  and  shallow  angle 
approximation 

Atmospheric  temperature  and  wind  are  mainly  functions 
of  height.  Therefore,  to  study  average  signal  amplitude  and 
long-period  signal  variations,  a  one-dimensional  model  can  be 
used.  The  signal  penetrating  height  depends  on  the  wind  shear 
and  the  temperature  gradient  in  the  atmospheric  boundary  layer 
as  well  as  the  source-receiver  distance.  With  a  horizontal 
travel  distance  of  9.2  km,  the  signal  penetrates  almost  to 
the  top  of  the  boundary  shear,  which  is  200  m  to  600  m  above 
the  surface.  The  signal  ray  angle,  c{  ,  which  the  ray  makes 
with  the  horizontal  surface,  is  always  small  (less  than  .3 
rad).  To  take  advantage  of  this  fact,  a  shallow  angle  ap¬ 
proximation  is  used  in  evaluating  ray  integrals. 

As  a  result  of  one-dimensional  model  and  shallow  angle 
approximation,  one  can  write 


cos  di  -  1  - 

0t2/  2, 

(la) 

COS  do~  1  ~ 

cVo2/2' 

(lb) 

1/C  a  1/CD  - 

<c-c0)/c02( 

(lc) 

by  Taylor's  expansion  and  retaining  terms  up  to  the  first 
order.  Then  the  substitution  of  Eqs.  (1)  into  Eq.  (10)  in 
the  previous  section  transforms  the  Snell's  law  into  the 
simplified  form 

CX  2  ■  0(o2  -  2  <  =  -  "o  +  Vy  -  Vyo,/C°' 
where  the  subscript  “o*  denotes  the  ground  value. 


(2a) 
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3y  defining  the  effective  wind  as 

ve  =  e  -  cQ  +  vy  -  vyo,  (2b) 

the  Snell's  law  Is  conveniently  written  as 

d2  ■  d02  -  2V°o.  (3) 

Then  the  integral  of  ray  trajectory  becomes 

y  =  f  (dz/oO 
*  o 

•  fo  {  dZ  '  ^  ±  (  °*°  '  2  Ve/C°)!i)}'  441 

where  o(  is  positive  in  the  first  half  trajectory  (from  the 
source  to  the  turning  point)  and  is  negative  in  the  second 
half  trajectory  (from  the  turning  point  to  the  landing  point). 
The  penetrating  height,  H,  is  the  value  of  z  which  makes  ck 
vanish. 

Let  the  receiver  location  be  y  =  Y  and  z  =  Z  =  0.  Then 
the  horizontal  travel  distance  is 

(* 

Y  -  2  (dz/00 

Jo 

=  2  /H(  dz/(  0lQ2  -  2  Vco)*)  ,  (5) 

where  H  -  the  signal  penetrating  height. 

The  ray  propagation  height  averaged  over  the  horizontal 
travel  distance  is  called  the  average  ray  height,  i.e.. 


Hav  =  (2/Y) 


rY/2  /H 

zdy  =  (2/Y)  /  (zdz/oO. 

3  Jo 
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3-3.  Fundamental  effective  wind  profiles  in  the  atmospheric 
boundary  layer 

The  diurnal  variations  of  the  average  signal  amplitude 

(see  Section  4-5  for  the  experimental  evidence)  suggest  6 

fundamental  effective  wind  profiles  which  are  applicable  in 

different  times  of  the  day.  They  ares 

1)  The  negative  shear  profile. 

Vs—  SZ,  (1) 

e 

where  S  is  a  positive  constant.  It  occurs  from  about  noon 
until  sunset.  No  signal  transmission  is  possible  with  this 


profile. 

2)  The  parabolic  profile. 

Ve  -  (2Vra  z/zm)  (  1  -  z/(2zm)] 

s  So  2  (1  -  z/(2zm)  J  , 


(2) 


where  zm  =  the  shear  vanishing  height, 

V_  =  the  wind  at  z_, 
m  mr 

SQ  =  the  surface  shear. 

It  is  applicable  mostly  at  night,  when  the  penetrating  height 
is  close  to  the  shear  vanishing  height  and  the  focusing  factor 
is  small.  The  parabolic  profile  can  be  derived  by  Taylor’s 
expansion  of  the  effective  wind  at  the  shear  vanishing  height, 


i.e.. 


V 

e 


(3) 


Since  the  effective  wind  vanishes  at  the  surface,  the  coef- 
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ficient  of  the  second  term  is  determined  as 


V 

e 

2 


m 


(4) 


which  is  a  negative  curvature.  Then  the  substitution  of  Eq. 
(4)  into  Eq.  (3)  gives  Eq.  (2). 

3)  The  logarithmic  profile. 

Ve  =  \  m  (z/za)  (5) 

where  zQ  =  roughness  length  =  1  m, 

-  the  wind  at  3  m- level. 

This  profile  is  predicted  assuming  a  strong  turbulent  boun¬ 
dary  shear  in  neutral  air  (Lumley  &  Panofsky  1964,  pp.103, 
Thuillier  &  Lappe  1964).  It  is  applicable  for  short-distance 
propagation  (Y  ^  2  km)  in  the  daytime  and  for  longer  distances 
at  night.  The  focusing  factor  is  usually  small. 

4)  The  linear  profile. 

Ve  =  Sz.  (6) 

It  is  applicable  at  night,  when  the  penetrating  height  is  in 
the  lower  part  of  a  strong  boundary  shear  and  the  focusing 
factor  is  about  1. 

5)  The  positive  curvature  profile.  This  has  three  types 
according  to  the  surface  shear,  which  can  be  positive,  zero, 
or  negative. 

Ve  =  SQ  z  [  1  +  z/(2zp)  |  ,  ^fype  a)  (7a) 

ve  =  so  z2/<2zp) 


# 


(Type  b) 


(7b) 
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Ve  =  3o  2  (  ~  1  +  z/<2zp)  )  ,  (Type  c)  (7c) 

where  SQ/Zp  is  the  curvature  and  the  surface  shears  are 
+  SQ,  0  and  -  S0.  These  profiles  are  applicable  when  the 
signal  penetrates  the  lower  part  of  strong  boundary  shears. 

The  focusing  factor  is  greater  than  1. 

6)  The  elevated  shear  profile.  This  is  a  two- layer 
model.  The  lower  layer  is  one  of  the  positive  curvature  pro¬ 
files,  while  the  upper  layer  3s  a  shear  with  negative  curva¬ 
ture.  The  interface  of  the  two  layers  is  an  inflection  point. 
The  elevated  shear  is  believed  to  be  responsible  for  the 
strong  signals  which  are  observed  in  the  beginning  as  well  as 
at  the  end  of  the  daily  receiving  period.  (Pig.  4-4.1). 

The  signal  properties  for  parabolic,  logarithmic,  linear, 
and  positive  curvature  profiles  can  be  evaluated  analytically 
(Appendix  2).  The  results  are  shown  below  with 
Y  =  the  travel  distance, 

H  =  the  penetrating  height, 

Hav  =  the  average  ray  height,  and 
f  =  the  focusing  factor. 

PARABOLIC  PROFILE 

Trajectory: 

2  —  ^  ^  0/o/  0(m)  sinh  y/zm)  +  1  —  cosh  (Q^y/z^)  ] 

-  =  (2V<V  tanh'1  <  <*o/o(m) 

H  =  zm  {  1  -  (  1  -  (  tfo/0(m)2  ]  h  } 


Hav  “  *m  (  1  *  <  <Vo/0(m)/tanh-1  (  o/o/ci  m)  ) 


«' ’m  *  <s0  V’o^  *  ‘2  V’c,)*5 


LOGARITHMIC  PROFILE 
Trajectory: 

y  =(1^  zQ  exp  (GQ2)/<a£  ]  (erf  (Gq)  ?  erf 

Y  =  2  IT3*  z0  exp  (Gq2)  erf  (GQ)/ 

K  =  zQ  exp  (Gq2) 

Hav  =  zc  exp  (Gc2)  erf  (2^  GQ)/  (2^  erf  (Gq) ] 

h  2  * 

IT  exp  <G0  )  erf  (GQ) 

.  2  G0(l  +  IT*  Gq  exp  (Gc2)  erf  (GQ)] 

Go  =  <V*/<2  V£  J35 

°(e  ®  (2  v£  /Co)35 

G  =  (go2  -  In  (z/zQ)  j35 

LINEAR  PROFILE 
Trajectory: 

<y  -  2o0fo/S)'!  =  (2  C./S)  [  -Q 0fo2/ ( 2S )  - 

V  =  2  ays 


Hav  =  (2/3)  H 
f  =  1 

positive  CURVATURE  TYPE  A  (positive  surface  shear) 
Trajectory: 

z  =  ZP  (  (0^o/^p)  sin  <dpy/V  ‘  1  *  GOS  (0^Py/zp) 
Y  =  (2  zp/ 0(p)  tan”1  (ol^Ol p) 

H  =  zp{(l  +  {Of^tfp)2)53  *  1  ) 

Kav  *  2p  (  <  O^o/OfpJ/tan*1  ( Ol^dl p)  -  1  ) 


f  - 

0(pY] 

2  z 

* 

cot 

°<PY] 

9  7 

+  tan 

CXpY 

9  fj* 

■ 

4 

► 

[  P  J 

£  z 

P  J 

^p  ~  <so  zp//Go^ 


POSITIVE  CURVATURE  TYPE  3  (zero  surface  shear) 
Trajectory: 

z  =  <#o  zp/0(p)  sin  (0(py/zp) 

Y  =  7T  zp/o/p 
^  —  &  o  zp/  0(  p 
Hav  =  ( 2/7T  )  H 
f  =  CO 

POSITIVE  CURVATURE  TY'PE  C  (negative  surface  shear) 


Trajectory: 
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z  =  zp  (  <  Ct0/Ctp)  sin  (  0(p  y/Zp)  +  1  -  cos  (  <tfpy/zp)  ) 

Y  =  (2  zp/o<p)  (  TT/2  +  cot"1  (Ota/Otp)  ] 

H  =  2P  {  1  +  (  1  +  i0(c/0(p)2}h} 

Hav  =  zp  {  1  +  ( 0(o/0(p)/  (  7T/2  +  cot”1  (  o(0/o(p)  ]  j 
f  =  (0^/0^)*  (l  +  (  cIq/  Ofp) 2  [  TT/2  ♦  cot"1  (<tf0/crp)]^ 

When  the  focusing  factor  is  small  at  night,  either 
parabolic  or  logarithmic  profile  is  possible.  The  best  way 
to  distinguish  between  these  two  would  be  to  have  several  re¬ 
ceivers  located  at  different  distances  and  aligned  with  one 
source,  and  to  measure  the  focusing  factor  as  a  function  of 
distance.  Fig.  3-3.1  compares  the  characteristics  of  para¬ 
bolic  and  logarithmic  profiles.  Both  profiles  have  a  wind 
of  10  m/s  at  400  m.  The  focusing  factor  of  the  logarithmic 
profile  drops  very  fast  within  the  first  kilometer,  and  then 
decreases  very  slowly,  having  a  value  near  .3  for  the  next 
8  km.  The  focusing  factor  of  the  parabolic  profile  is  almost 
unity  in  the  first  kilometer  (behaving  like  that  of  the 
linear  profile^)  and  then  decreases  fast  with  distance. Fig. 
3-3.2  depicts  f-contours  in  the  log  zm  vs.  log  diagram 
for  Y  =  9.2  km  and  0o  =  331  m/s,  which  helps  one  determine 
the  appropriate  parabolic  profile.  The  penetrating  heights 
and  average  ray  heights  as  fractions  of  zm  of  the  parabolic 
profile  are  plotted  vs.  f  in  Fig.  3-3.3. 

The  focusing  factor  of  the  linear  profile  is  almost 
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unity,  decreasing  slowly  with  the  travel  distance  (Appendix 
2).  The  focusing  factor  of  the  positive  curvature  profile 
is  greater  than  unity.  Travel  distances  of  positive  curva¬ 
ture  profiles  are  restricted  within  the  following  limits, 
which  are  also  caustic  distances  except  for  the  origins 

0<Y  <  TT  zp  / 0(p*  (Type  a) 

V  =  IT  zp  /0(pi  (Type  b) 

lTzp  /<*p<Y  <  2lTzp/°tp-  (Type  c) 

In  reality,  the  finite  thickness  of  the  atmospheric  layer, 
where  the  profile  is  valid,  limits  the  travel  distance. 

Fig.  3-3.4  shows  schematically  the  signal  penetrating 
heights  and  focusing  factors  vs.  the  travel  distance  for  the 
elevated  shears.  An  elevated  shear  predicts  a  large  focusing 
factor  when  the  penetrating  height  is  near  the  point  of  in¬ 
flection  of  the  profile.  Elevated  shear  type  b  or  c  predicts 
a  skip  distance  Ymini  which  is  also  the  caustic  distance. 

Ymini  can  be  estimated  by  assuming  the  elevated  shear  to  be 
constant  (Appendix  2).  Thus  one  gets  for  type  b, 

^mini  =  4  (~o  (^ 

where  =  the  thickness  of  the  lower  constant  velocity  layer, 

S2  -  the  constant  elevated  shear, 
and  for  type  c, 

vmini  =  2  12  =0  Vsl>%  (  U  +  Si/s2)2  -  1  J  \  UO) 


Fig.  3-3.4  Penetrating  heights  and  focusing 
factors  with  elevated  shears 


PROFILES 


PENETRATING  HEIGHTS 
vs.  TRAVEL  DISTANCES 


FOCUS IN i  FACTORS  vs. 
TRAVEL  DISTANCES 


Tvpe  b  or  c 
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3-4.  Amptitude  variations  and  doppler  shifts  due  to  parabolic 
profile  fluctuations 

A  linear  theory  is  developed  to  calculate  the  variations 
of  both  amplitude  and  frequency  due  to  fluctuations  of  the 
parabolic  effective  wind  profile.  With  this  theory,  the  ef¬ 
fects  of  individual  profile  parameters  on  the  signal  variation 
are  linearly  added.  The  parameters  for  parabolic  profile  are 
SQ  (surface  shear)  and  zm  (shear  vanishing  height).  The  re¬ 
sults  of  the  calculation  (Appendix  3)  are  the  following.  The 
predicted  time  rate  of  fractional  change  of  amplitude  due  to 
the  profile  fluctuation  is 

£a  ■  (St  (ln  V  )  <-  1  ♦  a-  <« 

where 

^  (In  SQ)  =  the  time  rate  of  fractional  change  of  SQ, 

Q  =  d  (In  zm)/d  (In  SQ) 

=  the  ratio  of  fractional  variations  of  z  and 

o  Itt 

bo' 

and 

A=  <’4){(soY/(2Co0(o))[l+  (  0/o/  0fm> 2  ]  -  1  j 

=  the  ratio  of  fractional  variations  of  signal  amplitude 

and  z  . 
m 

The  predicted  fractional  doppler  (shift)  due  to  the  profile 
fluctuation  is 

*d  '(St  <ln  so)  ]  <DSo  +  a  Dzm),  (2) 

where 
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ESo  =  (  \V2:.l)[(so  Y/(2C040)  )  (  1  ♦<  0to/dm)2}-  1  } 

=  the  doppler  due  to  variations  of  SQt 

Bzm  =  (  3  *m  ^£/‘2co>){(sov/<2co°(o)  H1  '<  <X</<*m)2'3  )  ‘  1  ] 

=  the  doppler  due  to  variations  of  zm. 

One  can  infer  the  profile  fluctuation  from  the  observed 
long-period  signal  variations.-  Let  the  observed  time  rate  of 
fractional  change  of  amplitude  due  to  the  profile  fluctuation 
be  F&,  the  associated  fractional  doppler  be  F d,  and  the  ratio 
of  Fa/Fd  be  F.  Then,  by  equating  the  observed  Fa  and  F^  to 
the  predicted  fa  and  f^  respectively,  one  obtains  the  inferred 
ratio  of  fractional  variations  of  zm  and  SQ  as 

Q  =  (FDSo  +  A)/(-  FDzm  *  A),  (3) 

and  the  inferred  time  rate  of  fractional  change  of  sQ  as 

It  <ln  so)  =  (Fa/A)/(-l*Q)  =  Fy(Dso  +  Q  Dm).  (4) 

Finally,  one  gets  the  inferred  profile  fluctuation  (peak  am¬ 
plitude)  as 

AVe  =  (r/4)  (  ft  (ln  SqH  jsoz(  X  -  (X  -  Q)z/(2  zm)]},  (5) 

where  'X  =  the  period  of  variation. 

Schematic  diagrams  in  Fig.  3-4.1  illustrate  all  pos¬ 
sible  variations  of  parabolic  profiles  according  to  the  linear 
theory.  The  eminent  properties  of  these  variations  are  listed 
in  Table  3-4.1,  where  zarn  =  2^/(1  -  Q)  denotes  the  height  of 
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Table  3-4.1  characteristics  of  parabolic  profile  fluctuations 


Fluc¬ 

tuation 

type 

0( 

P 

r 

i 

£ 

Q 

♦  ♦ 

♦  -  + 

+  ♦  ♦  ♦ 

+  4-  ♦  + 

+  -  +  + 

F 

zam 

zao 

(Ap)(A  V) 

Remarks: 

1.  The  fluctuation  type  is  determined  by  values  of 
Q  and  F  as  shown  in  Fig.  3-4.1. 

2.  Q  is  the  ratio  of  fractional  variations  of  zm 
and  SQ  of  the  profile,  i.e.,  d(ln  zm)/d(ln  SQ). 

3.  F  is  the  observed  ratio  of  the  time  rate  of 
fractional  amplitude  change  to  the  fractional 
doppler  shift. 

(1  -  Q)  is  the  height  where  the 
wind  is  relatively  maximum. 

5.  zao  =  2  zam  is  the  height  where  the  fluctuating 
wind  is  zero. 

6.  (Ap)(AV)  is  the  product  of  the  signal  amplitude 
change  and  the  wind  change  in  the  boundary  layer. 
When  alternative  signs  are  possible,  the  positive 
sign  is  for  z  >  zao  and  the  negative  sign  for 

z  <  2ao» 


fluctuating 
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a  relatively  maximum  fluctuation  and  zao  =  2  zam  denotes  the 
zero  fluctuation  height.  Also  listed  in  Table  3-4.1  is  the 
sign  of  the  product  of  the  signal  amplitude  change  and  the 
wind  change,  i.e.,  (&  p)  (Av)»  For  fluctuation  type  ck  ,  it 
is  positive,  wwi  h  means  that  an  amplitude  increase  accom¬ 
panies  wind  increases  at  all  levels  in  the  boundary  layer. 

For  types  ft  and  "Jf,  (A  p)  (AV)  is  negative.  For  types  $ 
and  £  ,  the  sign  depends  on  the  height  where  the  wind  is  re¬ 
corded.  The  observed  signs  of  (Ap)  (  A  V)  and  F  should  help 
identify  the  possible  types  of  profile  fluctuation  which  cause 
the  signal  variations. 

For  the  theory  of  profile  fluctuations  to  be  applicable, 
the  period  of  signal  variations  must  be  long  enough  to  insure 
a  practical  uniformity  of  associated  wind  variations  over  the 
source-receiver  line,  since  the  theory  is  based  on  a  one¬ 
dimensional  atmosphere  model.  If  the  fluctuating  wind  system 
drifts  with  the  average  wind,  which  blows  at  5  m/s  along  the 
9  km  source-receiver  line,  for  example,  then  the  period  of 
observed  profile  fluctuations  must  be  at  least  30  min.  If 
the  profile  fluctuation  wind  is  caused  by  gravity  waves,  which 
have  a  phase  velocity  component  of  50  m/s  along  the  9  km 
source-receiver  line,  for  example,  then  the  signal  period 
must  be  at  least  3  min.  The  periods  of  all  doppler  shift  data 
in  this  experiment  are  not  more  than  15  min.  Cne  might  hope 
to  get  profile  fluctuations  due  to  gravity  waves.  Correlations 
among  the  probe  signal  variation,  the  microbarograph  array 
data,  and  the  wind  records  were  searched  for,  but  no  definite 
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correlation  was  found.  A  theory  seems  to  be  needed  to  pre¬ 
dict  characteristics  of  gravity  waves  in  atmospheric  boundary 
shears,  and  is  not  yet  available. 

In  order  to  get  a  rough  idea  of  what  profile  fluctua¬ 
tions  might  possibly  be  inferred  from  the  field  data,  the 
long-period  signal  variations  on  Nov.  21-22,  1970  (Table 
A-1.6)  were  analyzed,  and  the  following  types  of  profile 
fluctuations  were  obtained: 

1)  Profile  fluctuation  type  S  (F  =  -  32  sec”^,  Q  =  -  1.4), 
for  which  the  maximum  wind  fluctuation  has  a  peak  ampli¬ 
tude  of  .32  m/s  at  123  m. 

2)  Profile  fluctuation  type  £  (F  =  +  32  sec~^,  3  =  -  2.9), 
which  has  the  maximum  wind  fluctuation  of  .12  m/s  (peak 
amplitude)  at  76  m. 

The  average  ray  height  that  night  was  168  m. 
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3-5.  Short-period  fluctuations 

When  fluctuation  periods  are  from  about  .5  to  8  min., 
the  turbulence  scale  is  small  compared  with  the  signal  travel 
distance,  and  three-dimensional  wind  fluctuation  models  should 
be  considered.  Although  there  has  been  much  work  in  the  li¬ 
terature  about  wave  propagation  in  a  random  medium  (e.g. 
Chernov  1960,  Tatarski  1961),  signal  fluctuations  in  the  ran¬ 
dom  medium  with  average  shear  structures  do  not  seem  to  have 
been  studied  before.  In  this  section,  amplitude  variations 
and  phase  shifts  of  the  signal  propagating  in  turbulent  boun¬ 
dary  shears  according  to  ray  theory  will  be  presented.  To 
apply  the  formulas,  one  should  beware  of  limitations  of  the 
theory.  Scattered  rays  due  to  atmospheric  inhomogeneities 
may  interfere  with  the  direct  ray  when  the  propagation  dis¬ 
tance  relative  to  inhomogeneity  size  scales  exceeds  a  cer¬ 
tain  limit.  After  analyzing  the  field  data  of  probe  measure¬ 
ments,  one  will  see  that  diffraction  phenomena  are  not  ne¬ 
gligible  in  this  experiment.  For  transmissions  to  the  9.2  km 
site,  the  wave  parameter,  which  serves  as  the  criterion  divi¬ 
ding  ray  and  diffraction  regimes  (Chernov  1960,  p.  74),  is 
found  to  be  of  the  order  of  10.  The  field  data  is  on  the 
diffraction  side,  not  far  from  the  border  of  ray  and  diffrac¬ 
tion  regimes.  However,  no  full-wave  solution,  which  covers 
both  ray  and  diffraction,  has  yet  been  available  for  signal 
propagation  in  a  turbulent  boundary  layer.  Therefore,  the 
field  data  of  short-period  fluctuations  will  be  interpreted 
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by  the  ray  theory  formulas,  and  possible  corrections  for  dif¬ 
fraction  will  be  estimated  (section  4-1). 

Details  of  deriving  ray  theory  formulas  are  described 
in  Appendix  4.  The  shallow  angle  approximation  is  used.  The 
atmosphere  is  assumed  to  have  a  one-dimensional  effective 
shear  and  three-dimensional  wind  fluctuations,  whose  correla¬ 
tion  function  has  a  three-dimensional  Gaussian  form,  i.e. 


where  a„,  a  and  a„  are  the  correlation  scales  alona  three 
X'  y  z 

coordinate  axes.  The  results  of  the  derivation  are  the 
following. 

The  mean  square  phase  shift  is 

(  A<}»2  «  <7 TH  CO2  ay  Y/c4)  (Ave)2  (2) 

where  =  the  angular  frequency  of  the  signal, 

a  =  the  scale  alona  the  source-receiver  line, 

y 

C  =  the  average  sound  speed, 

- - 

(Ave)  -  the  mean  square  effective  wind  fluctuation. 

With  the  shallow  angle  approximation,  phase  shifts  due  to 
effective  wind  fluctuations  do  not  depend  on  average  effective 
wind  profiles.  The  phase  shift  is  directly  proportional  to 


ci-tcCulVC 


W  JwilVA  Am  AUW  UUU  W-JuOli  UUU-  OUUu*-w  —  -  -  -  ... 


and,  is  independent  of  a_  and  av. 

4  A 
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The  mean  square  fractional  amplitude  variation  varies 
with  the  average  effective  wind  profile.  For  the  linear  pro¬ 
file, 

Tip)7  /  P2  =  (  THJ7  a/C2  ) . 

.  (  (7 Th  Y3/10)/az4  +  (7T35  Y3/15)/(az2  ax2) 

+  (7Th  y3/10)  /  ax4]  ,  (3) 

which,  for  equal  ax,  ay,  and  az,  becomes 

(Ap)2  /  P2  =  (4/15)  IT*  Y3  (  (Ave)5  /  a3  C  2  )  ,  (4) 

where  a  =  av  -  a„  -  a_.  In  comparison  with  a  plane  wave 
source  (Chernov  1960),  the  mean  square  amplitude  fluctuation 
of  a  point  source  in  a  linear  profile  is  10  times  smaller. 

The  result  of  Eq.  (4)  is  the  same  as  for  a  spherical  wave  in 
a  uniform  average  medium  (Tatarski  1961).  For  the  parabolic 
profile, 

(Ap)2  /  P2  =  (  <Aver  ay/C2  ]  . 

•  (  Kl/az4  +  K2/(az2  ax2)  +  Vax4  j  •  (5) 

%  3 

The  coefficient  K3  (  =7T  Y  /1C)  is  the  same  as  for  the  linear 
profile.  Coefficients  and  K ^  ate  essentially  functions 
of  the  focusing  factor  for  a  given  Y  (the  travel  distance). 
When  the  focusing  factor  becomes  unity,  and  become  the 
corresponding  coefficients  of  the  linear  profile,  i.e., 

Kx  Y3/10, 


m  i".  •  w  * 


ji.JUMJ,  1 1 1)\{§XM\ |» !  WfPWWOT 


64 


and 


TT^  Y3/15,  as  f  - >  X. 


Fig.  3-5.1  plots  the  coefficients  K^,  K3  as  well  as  the 
sum  of  the  three,  K„ ,  vs.  the  focusing  factor  f. 

•"3 

Kow  the  way  of  inferring  the  magnitudes  and  scales  of 
effective  wind  fluctuations  from  the  field  data  is  described. 
The  turbulence  scale  or  average  eddy  size  is  given  by  the 
integral  (Taylor  1935) 


L  =  J  f 

Jo 


R  (s)  ds, 


(6) 


where  L  =  the  integral  scale, 

.i(s)  =  v(r)v(r+s)  /  v2 

=  the  normalized  one-dimensional  correlation  function 
of  wind  fluctuations, 

s  =  the  separation  coordinate. 

For  the  three-dimensional  Gaussian  correlation  function,  the 
turbulence  scales  are 


=  nr  /2)  ax 


b  =  (7TV2)  a 

y  y 


Lz  =  (lTh/2)  az 


(The  horizontal  scale  perpendicular 
to  the  source-receiver  line),  (7a) 


(the  horizontal  scale  along  the 
source-receiver  line),  (7b) 


(the  vertical  scale). 


(7c) 


From  the  field  data,  one  obtains: 

Feriod  of  signal  variation  =?  (sec) 
Fractional  rms  amplitude  variation  =  (T^/p 
rms  phase  shift  =C^(rad) 


3 

J 
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Focusing  factor  =  f 

some  data  gi  'es  the  doppler  instead  of  the  phase  shift 
.  The  conversion  formula  is 

a4>  =  T  A  JS  .  (?) 

An  appropriate  parabolic  profile  can  be  chosen  from 
the  focusing  factor  and  the  synoptic  weather  information. 

1o  estimate  the  horizontal  turbulence  scale,  one  assumes 
that  the  turbulence  is  horizontally  isotropic.  In  other  words, 
the  horizontal  turbulence  scale  along  the  source-receiver  line, 
i.e.,  by,  is  assumed  to  be  equal  to  the  horizontal  turbulence 
scale  perpendicular  to  the  source-receiver  line,  i.e.,  Lx. 

•since  short-period  signal  variations  are  mainly  caused  by  the 
turbulent  air  motion,  drifting  with  the  average  wind  across 
the  source-receiver  line,  the  horizontal  turbulence  scale  is 
approximately  the  product  of  the  period  of  signal  variation 
and  the  cross  wind  component,  i.e., 


*  ax  =  (vrr')T:  Vi;av 


(9) 


where  V  =  the  magnitude  of  the  horizontal  wind  comconent 
xav 

perpendicular  to  the  source-receiver  line  at  the  average  ray 
heig.-  (nav). 

Then  the  rms  wind  fluctuations  is 


/  (Tl'-CUa^  V‘J) 


(10) 


jy  taking  the  ratio  of  phase  and  amplitude  variations 
given  by  the  field  data,  one  nets  from  ..nr,.  ( 2 )  and  ( 5 } 
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o$/  ( oyp)  =  or*'  ^ 


7  (Va24  +  V(az2  ax2)+K3/ax4]  ^  (11) 


from  this  equation,  one  can  determine  a7  and  ax,  if  the  ratio 


of  a7a,  (=M)  is  _;iven,  i.e., 

/T  z 

a2  “  3  ^  0$  V(lT1/s  Y*  (  <J^/F)%]  (12a) 

(12b) 


a  =  i-i  a„, 
x  z 


where  K.,..  =  K,  +  K0  /  M2  +  K,  /  K** . 
Si'';  1  Z  3 


»»4 


(12c) 


First  of  all,  isotropic  eddies  (av  =  a,r  =  a  )  are  assumed. 

Cne  puts  M  =  1,  and  evaluates  ax  (=a2)  by  Eqs.  (12).  The 

ax  thus  obtained  is  supposed  to  be  of  the  same  order  as  a^, 

which  is  determined  by  Zq.  (9).  otherwise,  horizontally 

isotropic  eddies  (a  =  a  ?£  a  )  are  assumed.  Cne  should  try 

x  y  z 

various  values  of  \u  until  a,  is  of  the  same  order  as  a,,. 

a  y 

From  the  experimental  results,  the  values  of  a_  are  found  to 

Cu 

be  either  of  the  same  order  of  or  smaller  than  the  values  of 
a  ,  i.e.,  ^1  (Section  4-3).  As  i-i  inceases,  the  contri- 

X 

butions  of  I'_  and  K,  to  decrease  very  fast,  and  that  of 
remains  the  same.  -Vhen  M  2>  3,  almost  equals  F^. 


remains  the  same 
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CHAPTER  4.  EXPERIMENTAL  RESULTS 

Most  of  the  field  data  were  received  at  9.2  km  to  the 
SE  with  the  Daniels  pipe  as  the  noise  reducing  device.  For 
longer  travel  distances/  signal  has  been  received  at  11  km  to 

the  NE  and  at  20  km  to  the  E  on  exploratory  field  trips  using 

wind  screens,  but  the  information  about  signal  variabilities 
is  not  useful  because  of  poor  signal-to-noise  ratio.  For  a 
travel  distance  of  4.5  km  to  the  E,  there  is  some  good  field 
data  obtained  with  only  wind  screens. 

There  is  a  small  hill  between  the  source  and  the  re¬ 
ceiver  at  9.2  km  to  the  SE  as  shown  in  Fig.  4.1.  The  direct 

ray  must  have  an  initial  elevation  angle  of  at  least  .37  deg 

in  order  to  be  received.  Transmission  experiments  to  the  9.2 
km  receiver  were  frequently  carried  out  from  June  29  through 
Dec.  30,  1970.  The  field  data  collected  during  this  half  year 
has  two  types  of  signal  variation.  The  variations  with  shorter 
periods  of  .5  to  8  min  have  both  amplitude  and  doppler  (or 
phase)  information.  The  variations  with  longer  periods  of 
.5  to  6  hours  have  only  amplitude  information. 

Section  4-1  gives  examples  showing  how  to  analyze  short- 
period  field  data  to  infer  wind  fluctuations  and  turbulence 
scales,  and  estimates  possible  errors  induced  by  the  inac¬ 
curacy  of  the  ray  theory  and  analyzing  methods.  In  order  to 
investigate  diffraction  effects,  working  formulas  are  derived 
for  predicting  amplitude  variations  and  phase  shifts  of  the 
signal  propagating  in  turbulent  boundary  shears  with  wave 


Fig.  4.1  Map  of  the  experiment  area 


°^50  oR4.5 
o«w3  (-12m) 
6m) 


(Om) 


S:  signal  source 

W3 :  3m- level  wind  sensor 

w50s  50m- level  wind  sensor 

R4  5:  receiver  4.5  km  to  the  E  (or  84°) 

Rg  2s  receiver  9.2  Jem  to  the  SE  (or  134°) 
(*m)j  *  is  the  elevation  above  Rq  in  meters 
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parameters  of  the  order  of  10.  Sections  4-2  and  4-3  show  the 
characteristics  of  wind  fluctuations  and  the  shape  of  eddies 
in  the  atmospheric  boundary  layer.  Section  4-4  discusses 
possible  signal  variations  due  to  gravity  waves.  The  most 
notable  long-period  phenomenon  of  the  atmospheric  boundary 
layer  is  the  diurnal  variation.  This  is  presented  in  Section 
4-5,  which  is  also  the  experimental  evidence  for  fundamental 
effective  wind  profiles  discussed  in  the  previous  chapter. 

The  wind  speeds  at  3  m-  and  50  m- levels  were  usually 
recorded  when  the  probe  was  operated.  The  locations  of  the 
anemometers  are  all  near  the  signal  source  as  shown  in  Fig. 
4-1.  The  wind  direction  and  surface  temperature  can  be  es¬ 
timated  from  the  hourly  weather  reports  of  Bedford,  Mass., 
which  is  about  20  km  to  the  SE.  Regular  upper  air  soundings 
are  taken  twice  a  day  (0700  and  1900  EST)  at  Albany,  N.Y., 
Nantucket,  Mass,  and  Portland,  Maine,  which  sometimes  indicate 
the  air  mass  characteristics  overhead.  These  weather  stations 
are  about  180  km  away  from  the  experiment  area.  The  wind  pro¬ 
file  in  the  lowest  1  km  layer  can  sometimes  be  obtained  from 
pibal  soundings  at  Boston,  Mass.,  which  is  about  40  km  to  the 
SE. 
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4-1.  Short-period  amplitude  variations  and  doppler  shifts 
(or  phase  shifts)  of  the  received  signal  and  their 
analysis 

The  phase  shifts  and  their  associated  amplitude  varia¬ 
tions  of  the  probing  signal  with  periods  of  .5  to  8  min. 
allow  one  to  infer  effective  wind  fluctuations  and  turbulence 
scales.  The  effective  wind,  as  a  convenient  concept  for  sig¬ 
nal  propagation  in  the  atmospheric  boundary  layer,  simplifies 
the  theory.  However,  with  a  single  receiver  as  used  in  this 
experiment,  one  can  not  distinguish  between  temperature  and 
wind  fluctuations.  To  the  first  approximation,  the  experi¬ 
mental  results  are  analyzed  by  assuming  that  signal  fluctua¬ 
tions  are  solely  due  to  horizontal  wind  fluctuations.  This 
assumption  is  based  on  the  following  three  reasons.  Firstly, 
the  temperature  fluctuation  is  less  efficient  in  affecting 
sound  propagation  than  the  wind  fluctuation.  One  degree  cen¬ 
tigrade  of  temperature  difference  corresponds  to  only  .6  m/s 
of  sound  speed  difference.  Secondly,  the  temperature  inho¬ 
mogeneity  is  very  efficient  in  producing  the  wind.  According 
to  the  thermal  wind  formula  (Haurwitz  1941,  p.  149),  for 
example,  one  degree  centigrade  of  temperature  difference  over 
one  km  of  horizontal  distance  gives  a  vertical  wind  shear  of 
.3  sec-^,  i.e.,  30  m/s  of  horizontal  wind  variation  per  100  m 
of  vertical  distance.  Therefore,  temperature  fluctuations 
seem  to  affect  the  signal  propagation  rather  indirectly  by 
producing  wind  fluctuatio'  s  than  directly  by  presenting  sound 
speed  fluctuations.  Thirdly,  one  will  see,  after  analyzing 
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the  field  data,  that  there  is  no  significant  difference  of 
the  wind  fluctuation  magnitude  between  the  night  data  and 
the  data  collected  in  early  stages  of  convective  periods 
(Table  4-1.1),  although  one  would  expect  more  temperature 
fluctuations  during  convective  periods. 

Three  examples  of  short-period  signal  variations,  one 
for  each  stage  of  equipment  condition,  will  be  presented  in 
this  section  to  show  how  wind  fluctuations  and  turbulence 
scales  are  inferred  from  the  field  data.  The  three  stages  of 
equipment  are,  as  mentioned  in  Section  2-2,  characterized  by: 
1)  a  4  cps-bandwidth  preamplifier  filter  and  recording  doppler 
(shifts),  2)  a  1  cps-bandwidth  filter  and  still  recording  dop¬ 
pler,  and  3)  a  1  cps-bandwidth  filter  and  recording  phase 
shifts. 

The  first  example  typical  for  the  first  stage  of  equip¬ 
ment  is  taken  from  the  June  30  data  (Fig.  4-1.1).  It  was  a 
turbulent  night  as  indicated  by  the  3  m- level  wind.  The  re¬ 
corder  for  50  m-level  wind  had  not  been  installed  at  that 
time.  There  was  a  thunder  shower  later  that  night  at  about 
0420  £ST.  The  signal  variations  shown  on  Fig.  4-1.1  must 
have  been  caused  by  the  boundary  layer  turbulence  before  a 
thunder  shower.  The  wind  was  blowing  from  235°,  and  makes  an 
angle  of  79°  with  the  source-receiver  line.  The  signal  source 
was  scheduled  to  be  on  for  14  min  and  off  for  2  min,  with  a 
cycle  every  16  min.  The  broken  line  on  the  signal  records 
shows  the  on  or  off  switching  time.  One  can  easily  read  the 
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signal-to-noise  ratio  by  comparing  AC  levels  of  on-  and 
off-periods.  The  average  signal-to-noise  ratio  is  about  2, 
which  is  enough  for  reliable  doppler  shifts.  The  average 
period,  amplitude  variation  and  phase  shift  are: 

=  52  sec,  (the  period  of  variation) 
g^/p  =  .147,  (the  rms  fractional  amplitude  variation) 

GJ  =  .312  rad.  (the  rms  phase  shift) 

The  parabolic  wind  profile  appropriate  for  that  night  has  the 
parameters : 

zra  =  300  (the  shear  vanishing  height) 

Vm  =  2.25  m/s,  (the  effective  wind  component  at  zm) 

H  =  201  m,  (the  signal  penetrating  height) 


Vxav  =  9'3 


(the  cross  wind  component  at  the 
average  ray  height) 


Kg  s  1.93  x  I0lj-  ra3. 


(the  coefficient  of  amplitude 
variation) 


Therefore  the  wind  fluctuation  and  turbulence  scales  are 
(j^  =  .138  m/s  (the  rms  wind  fluctuation) 


Lx  =  941  m. 


(the  horizontal  turbulence  scale  perpen¬ 
dicular  to  the  source-receiver  line) 


L  =  95  m. 


(the  vertical  turbulence  scale) 

(Wind  is  from  235°.) 

The  second  example  typical  for  the  second  stage  of  equip¬ 
ment  is  taken  from  the  Nov.  23  data  (Fig.  4-1.2).  It  is  also 
an  example  of  convective  data.  Although  a  passing  cold  weather 
front  obscured  some  of  the  regular  convective  phenomena,  the 
strong  signal  amplitude  between  0912  and  0934  EST  remains  to 
be  the  familiar  feature  associated  with  an  elevated  shear  of 
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the  early  convective  period  {Section  4-5).  The  doppler  shifts 
are  reliable,  because  the  variations  of  amplitude,  frequency 
and  AC  level  are  all  well  correlated  and  the  average  signal- 
to-noise  ratio  is  about  3,  which  is  greater  than  the  minimum 
requirement  of  1.  The  average  period,  amplitude  variation 
and  phase  shift  are: 

=  43  sec, 

0^/p  =  .050, 

GT£j>  =  .115  rad. 

The  appropriate  parabolic  profile  has  the  parameters: 

zm  =  360 

Vm  =  4.50  m/s, 

H  =  275  m, 

V  =  3.55  m/s, 
xav  ' 

Ks  =  1.77  x  10 11  m3. 

Therefore  the  wind  fluctuation  and  turbulence  scales  are: 

0^  =  .096  m/s, 

—  198  m, 

Lz  =  154  m. 

(Wind  is  from  250°. ) 

The  third  example  typical  for  the  third  stage  of  equip¬ 
ment  is  taken  from  the  Dec.  14  data  (Fig.  4-1.3).  Doppler 
shifts  were  integrated  and  recorded  directly  as  phase  shifts. 
The  phase  shifts  are  reliable,  because  signal-to-noise  ratios 
are  greater  than  1.  This  example  shows  large  signal  variations 
but  very  little  low-level  wind  fluctuations  of  comparable 
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periods.  The  signal  variations  seem  to  be  associated  with 
an  elevated  shear  in  the  winter.  No  significant  temperature 
inversion  in  the  lowest  1  km-layer  is  indicated  by  the  0700 
EST  upper  air  soundings.  The  average  period,  amplitude  varia¬ 
tion  and  phase  shift  ares 
T  =  72  sec, 

OJ/P  =  -223, 

<3^  =  .540  rad. 

The  suitable  parabolic  profile  has  the  parameters: 

zm  =  375  m, 
m 

vm  =  5.75  m/s, 

H  =  300  m, 

V  =  2.0  m/s, 
xav 

Kg  =  1.77  x  1011  m3. 

Therefore  the  wind  fluctuation  and  turbulence  scales  ares 


(J£  =  .463  m/s 
Lx  =  140  m. 


L  =  156  m. 
2 


(Wind  is  from  300  .) 

Finally  the  average  periods,  wind  fluctuations,  and 
turbulence  scales  derived  from  all  available  field  data  are 
summarized  in  Table  4-1.1.  The  details  are  listed  in  Ap¬ 
pendix  1. 

The  average  vertical  turbulence  scale  is  99  m  (Table 
4-1.1).  A  test  of  the  relative  importance  of  diffraction 
phenomena  is  the  wave  parameter,  which  is  defined  as  the  ratio 
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Table  4-1.1  Juramary  of  periods,  magnitudes,  and  scales  of 
horizontal  wind  fluctuations 


Date 

Time 

Period 

Magnitude 

Jcale 

of  the 

r 

bx 

(1970) 

day 

(sec) 

(m/s) 

(m) 

(m) 

5/8 

Kight 

188 

6/30 

Uight 

61 

.135 

982 

89 

8/15 

Is  ight 

111 

.143 

830 

53 

8/16- 

bight 

35 

.047 

402 

64 

17 

Jonv. 

29 

.107 

236 

78 

11/21- 

Kioht 

49 

.128 

219 

122 

22 

Jonv. 

49 

.057 

208 

106 

1.  ight 

455 

.229 

1364 

97 

11/23 

Jonv. 

43 

.090 

204 

153 

12/14 

i;  ight 

99 

Jonv. 

112 

12/19 

night 

104 

.181 

581 

126 

Jonv. 

105 

.191 

540 

104 

Kicht 

138 

.144 

730 

92 

Averaoe 

Jonv. 

68. 

.111 

297 

110 

Total 

Ill 

.131 

557 

99.2 

+.055 

+28 

Remarks : 


a)  The  travel  distance  is  4300  m  for  Kay  6,  1970  data, 
and  is  9200  m  for  the  remaining  data. 

b)  "Jonv."  for  the  time  of  the  day  refers  to  early  staces 
of  the  daily  convective  period  (approximately  from  1 
hr  after  sunrise  until  2-4  hrs  before  noon). 

c)  because  of  special  circumstances,  the  data  of  Kay  8 
and  December  14,  1970  does  not  give  meaningful  results 
about  wind  fluctuations  and  turbulence  scales  ('iables 
ii-1.1  and  *i.-l.  £ )  . 
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of  the  cross-sectional  area  of  the  first  Fresnel  zone  to  the 
area  scale  of  inhomogeneities,  i.e., 

D  =  2  /  (TT  a22) 

=  X  Y  /  (2  Lz2). 

The  wave  parameter  for  transmissions  to  the  9.2  km  is  found 
to  be  12.  Therefore,  the  diffraction  phenomena  are  not  ne¬ 
gligible.  cne  should  investigate  how  much  the  scattered  rays 
have  interfered  with  the  direct  ray. 

To  estimate  the  effect  of  diffraction  phenomena  on 
signal  fluctuations  in  the  boundary  turbulence,  one  can  modify 
the  formulas  of  intermediate  values  of  the  wave  parameter  which 
have  been  derived  for  a  plane  wave  (cbuchow  1953,  p.  167? 
Chernov  1960,  p.  83)  as 

(AfJ»2  =  E  {  1  +  (tan”1  C)  /  D)  (la) 

{  A  p)2  /  p2  =  E  (  1  -  (tan”1  r)  /  D J  (lb) 

where 

E  =  (  A  ve)a  TT95  ay  O)2  Y  /  (2  34), 

and  the  correlation  scales  perpendicular  to  the  source-re¬ 
ceiver  line,  az  and  ax,  are  here  assumed  to  be  equal.  In  the 
limit  of  ray  approximation,  i.e.,  for  small  wave  parameters, 
Eqs.  (1)  and  (2)  approach  Eqs.  (2)  and  (3)  in  Section  3-5  res¬ 
pectively  except  for  a  factor  in  the  formulas  of  amplitude 
variations,  which  varies  with  the  average  effective  shear. 

Thus,  in  the  ray  theory,  one  has 


T’r'} 


■w 
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(  =  2E  (2a) 

(^  p)2  /  p2  =  E  (KD)2  /  3,  (2b) 

where  D  is  the  wave  parameter  defined  above,  and  K  the  factor 
determined  by  average  shear  structures.  For  the  plane  wave, 

K  =  1. 

For  the  point  source  in  a  linear  effective  wind  profile, 

K  =  1  /  (10)^. 

For  the  point  source  in  a  parabolic  effective  wind  profile, 

K  =  (3KS  /  (8  TT3*  v3)  ) 

where  the  coefficient  Kg  for  Y  =  9200  m  can  be  read  from  Fig. 
3-5.1.  From  Eqs  (1)  and  (2),  one  can  have  generalized  for¬ 
mulas  of  intermediate  values  of  the  wave  parameter  for  signal 
propagation  in  the  turbulent  boundary  layer  as 

(p)2  =  E  (  1  +  (tan"1  D)  /  E)  (3a) 

( A  p)2  /  p2  =  E  {  1  -  ( tan”1  (KD)  )  /  (KD)  J .  (3b) 

These  formulas  should  allow  one  to  infer  more  nearly  correct 
values  of  wind  fluctuations  and  turbulence  scales.  From  Eq. 
(3a),  it  was  found  that  the  wind  fluctuation  has  been  underes¬ 
timated  by  using  the  ray  theory.  The  correction  factor  is 

C*v  ={2  /  (  1  +  (tan”1  D/D))}'5 


=  1.33  . 
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Therefore  the  wind  fluctuation  is  possibly  underestimated  with 
the  ray  theory  by  almost  33%.  .1  more  nearly  accurate  evalua¬ 

tion  of  turbulence  scales  calls  for  solving  the  equation,  from 
dqs .  (3), 

(.1,  -  1)  n  =  (R../IC)  tan"1  (XT)  +  tan"1  Z , 

where 

<6  -  2  /  (  0“ p2/?2)  (given  by  field  data), 

and  -v  is  given  by  the  average  boundary  shear  model.  The  ac¬ 
curacy  of  inferring  turbulence  scales  by  the  ray  theory  de¬ 
pends  on  R,.t  which  ranges  from  1.1  to  15.2,  averaging  at  5.3 
(Tables  A-l.l  to  .-i-1.9).  Ihe  correction  factor  for  diffrac¬ 
tion  is  found  to  be 

Ct  —  .34,  for  _  1  ■) 

:Lz  =  .98,  for  a.  =  5.3, 
and  :Lz  =  1*06'  for  ?'s  =  15»2» 

Therefore  the  error  of  inferring  turbulence  scales  by  the  ray 
theorv  is  crenerallv  within  5S  except  for  values  of  X  verv 

'*  J 

close  to  unity,  when  the  vertical  turbulence  scale  may  be 
overestimated  by  70b. 

The  results  of  data  analysis  can  have  errors  due  to  the 

inaccuracy  of  V..  (averace  horizontal  wine  component  oer- 
xav  -  - 

pendicular  to  the  source-receiver  line),  the  error  of  which, 
is  estimated  to  be  within  about  30  i.  Therefore,  the  errors 
of  inferred  wind  fluctuations  and  horizontal  turbulence  scales 
are  estimated  to  be  within  15  j  and  30’.  respectively,  .he  error 


of  inferred  vertical  turbulence  scales  caused  bv  the  inac- 


4-2.  The  inferred  horizontal  wind  fluctuations 

The  wind  fluctuations  inferred  from  the  probe  data  ap¬ 
proximately  represent  the  horizontal  wind  fluctuations  at  the 
average  ray  heights ,  which,  in  this  experiment,  range  from 
141  m  to  258  m,  averaging  at  192  m  (Table  A-1.10).  All  the 
field  data  were  collected  with  strong  boundary  shears.  Unless 
the  shear  is  elevated,  the  fluctuation  of  the  50  m- level  wind 
is  usually  a  good  indicator  of  the  boundary  turbulence.  The 
Nov.  21-22,  1970  data  is  a  typical  example  (Fig.  4-2.1).  The 
night  started  with  a  very  turbulent  boundary  layer  as  indicat¬ 
ed  by  the  wind  record.  It  became  generally  quieter  later  on, 
except  a  few  brief  periods  with  strong  turbulence  (e.g., 
0420-0510  EST).  The  average  ray  height  that  night  was  esti¬ 
mated  to  be  168  m.  The  inferred  rms  wind  fluctuations  (crv) 
vs.  the  rms  wind  fluctuations  monitored  at  50  m  (  (J^so)  are 
listed  in  Tables  A-l.ll  and  A-1.12  (Appendix  1)  and  plotted 
in  Fig.  4-2.2.  The  data  spreads  over  a  wide  range.  (Jv50  is 
about  1  to  14  times  (Tv.  If  a  linear  proportionality  exists 
between  Ov50  and  Ov*  the  50  m- level  wind  fluctuation  is 
roughly  3  times  the  inferred  wind  fluctuation  as  indicated  by 
the  eye-drawn  average  slope  in  Fig.  4-2.2.  Despite  the  wide 
spread  of  slopes,  the  inferred  wind  fluctuation  at  168  m  is 
clearly  smaller  than  the  wind  monitored  at  50  m. 

The  fact  that  the  wind  fluctuations  at  upper  levels  are 
smaller  than  those  at  lower  levels  in  the  boundary  layer  has 
been  demonstrated  by  direct  measurements  with  aircraft  (Bunker 


Fig.  4-2.2  rms  fluctuations  of  inferred  wind 
vs.  50  m-level  wind  on  Nov.  21-22 
1970 


Oi  short-period  data  (Table  A-l.ll) 
long-period  data  (Table  A-1.12) 


Wind  fluctuations  A 
inferred  from 
propagation 
fluctuations 
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1955,  Lettau  and  Davidson  1957),  and  with  anemometers  mounted 
on  high  towers  (Bysova  et  al  1965).  The  aircraft  measurements 
were  mostly  conducted  during  periods  of  strong  convection. 

The  horiz* tal  wind  fluctuations  detected  with  aircraft  ever 
the  North  Atlantic  Ocean  and.  over  the  C*Neill,  Neb.  plains 
all  roughly  -how  that  the  fluctuation  at  168  m  is  smaller  than 
that  at  50  m  by  a  factor  of  about  2.  The  Obninsk  300-m  tower 
data  includes,  a  stable  case  which  also  indicates  that  the  50 
m-f  lu>:  :.u<?  ...ion  is  greater  than  the  168  m-f  luctuation  by  a 
facte;  'v'  about  2.  7n  comparison  with  the  aircraft  and  tower 
measurements ,  the  inferred  wind  fluctuations  at  168  m  on  Nov. 
21-22,  1970  seem  to  be  slightly  smaller  but  of  the  same  order. 

The  rms  wind  fluctuation  averaged  over  all  the  available 
data  (Table  4-1.1)  is  .13  +  .06  m/s.  This  represents  the  hori¬ 
zontal  wind  fluctuation  at  about  192  m,  which  is  the  average 
Hav  as  listed  in  Table  A-1.10.  In  comparison  with  the  Obnink 
tower  data,  the  wind  fluctuation  inferred  from  propagation 
fluctuations  is  also  slightly  smaller  but  of  the  same  order. 

The  Obninsk  data  (Byzova  et  al  1965,  p.  79)  gi\  s  the  ratio  of 
^v/v300  at  192  m  as  about  .03  ♦  .006.  The  average  wind  at 
300  m  (v3qq)  for  this  experiment  is  estimated  to  be  about  10 
m/s.  Therefore  the  rms  wind  fluetation  at  192  m  should  be 
.30  +  .06  m/s.  Even  if  one  corrects  for  the  possible  underes¬ 
timation  by  the  ray  theory  (Section  4-1),  the  inferred  wind 
fluctuation  is  still  smaller  than  the  wind  fluctuation  measured 
with  meteorological  towers.  The  reason  for  a  smaller  wind 
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fluctuation  inferred  from  propagation  fluctuations  seems  to 
be  that  the  probing  signal  mainly  propagates  near  the  shear 
vanishing  height.  The  wind  near  shear  vanishing  heights  is 
expected  to  fluctuate  less  than  at  lower  levels. 
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4-3.  '.'he  inferred  eddy  sizes 

The  inferred  turbulence  has  horizontal  scales  of  204  m 
to  13C4  which  are  dictated  by  characteristics  of  the  re¬ 
cording  system,  the  inferred  vertical  turbulence  scales  range 
from  53  m  to  153  m  with  an  average  of  99  +  28  m.  therefore, 
the  eddies  seem  to  be  mostly  elongated  in  the  horizontal 
direction.  The  ratios  of  horizontal  to  vertical  scales  (L^/ 

>->')  are  slotted  vs.  horizontal  scales  doc)  in  Fig.  4-3.1.  The 
isotropic  eddies  have  scales  of  about  100  m.  As  the  horizontal 
scale  becomes  greater,  the  eddy  becomes  more  and  more  elongated 
horizontally. 

In  the  literature,  the  only  report  about  turbulence 
scales,  which  can  be  compared  with  the  results  of  this  ex¬ 
periment,  was  given  by  Lumley  and  Panofsky  (1964,  p.  196). 

They  analysed  the  data,  which  Singer  obtained  at  I-.rookhaven , 
and  concluded  that,  in  the  stable  air,  large  eddies  are  about 
4  times  as  long  (2000  m)  as  they  are  high  (500  m) ,  and  small 
eddies  are  more  nearly  isotropic.  They  did  not  mention  how 
small  their  isotropic  eddies  are.  However,  their  results 
agree  with  those  of  this  experiment  on  the  general  trend  that 
small  eddies  are  more  nearly  isotropic.  The  eddies  in  the 
boundary  shear  layer  as  inferred  by  the  probe  become  aniso¬ 
tropic  at  smaller  scales  than  what  Oinger  observed  at  Jrook- 
haven.  Jith  the  same  ratio  (L^/Lz)  of  4,  for  example,  the 
eddies  of  this  experiment  have  horizontal  scales  of  about 
42C  m  instead  of  2000  m  as  reported  by  bum ley  and  Panofsky. 


Reproduced  from 
best  available  copy. 


91 


In  other  words,  the  ecVj.ec  in  t>;o  ho: me' ary  shear  layer  nee:1, 
to  he  more  elongate''  horizontally,  the  share  or  the  eddies 
inferred  b**  the  orcbe  measurements  seems  to  have  been  modi¬ 
fied  hy  the  boundary  shear  in  such  a  way  that  they  arc  more 
elongated  horizontally  than  those  observed  hy  direct  measure¬ 
ments.  I  anofnky  *  linger  (1905)  and  tielke  y  Pan of sky  (1970), 
using  c'ata  from  various  locations,  calculated  the  phase  la' 
between  wind  fluctuations  at  different  heights.  ‘..he  upper 
portions  of  the  eddies  were  found  to  l^e  downwind  of  the  lower 
portions,  and  had  slopes  which  are  of  order  unity,  becoming 
more  nearly  vertical  with  height,  the  phase  lag  of  wind 
fluctuations  could  contribute  to  the  flattening  of  eddies, 
the  decrease  of  phase  lag  with  height  seems  to  imply  that  the 
wind  shear  '•'e creases  with  height  in  the  boundary  layer. 

It  is  also  interesting  to  investigate  the  relations 

among  the  vertical  turbulence  scale,  the  average  ray  height, 

anc  the  boundary  shear  thickness  ( zm) .  Fi<=  4-3.2  clots  the 

avera' e  ray  height  (riav)  vs.  the  vortical  scale  (L7).  r'he 

ratio  of  ij  /. H  is  found  to  be  .57  +  ,22.  Similar  1'’  one  gets 

z  av 

the  ratio  of  ij„/z,_  to  be  .29  +  .13.  In  other  words,  the  ver- 
tical  scale  of  the  eddies  is  roughly  57  i  of  the  average  ray 
height  and  29-i  of  the  7ooundary  shear  thickness. 


Fig.  4-3.2  H  vs.  L 

cl  V  ^ 
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4-4.  About  detecting  gravity  waves 

During  the  probe  measurements  over  the  9.2  km  travel 
distance  from  June  29  through  Dec.  30,  1970,  a  microbarograph 
array  was  continuously  monitoring  the  gravity  wave  pressure 
in  the  experiment  area.  The  microbarograph  records  were  com¬ 
pared  with  the  probe  signal  records.  No  significant  correla¬ 
tion  among  them  was  found.  According  to  a  simple  gravity 
wave  theory  (Madden  and  Claerbout  1968),  a  few  gravity  waves 
recorded  by  the  microbarographs  should  have  caused  detectable 
doppler  shifts  of  the  received  signal,  but  no  corresponding 
signal  frequency  variations  were  obtained.  The  actual  reason 
for  this  is  still  not  clear.  The  doppler  shifts  due  to  gravity 
waves  could  have  been  obscured  by  stronger  doppler  shifts  with 
shorter  periods  due  to  the  boundary  layer  turbulence.  Later 
in  the  experiment,  more  electronic  filtering  was  introduced 
to  the  receiving  system  in  order  to  show  long-period  signal 
variations  more  clearly.  Also  the  doppler  shifts  were  inte¬ 
grated  through  a  high-pass  filter  to  give  signal  phase  shifts. 
But,  unfortunately,  the  signal  transmitter  was  still  under 
the  constant-frequency  control  mode.  The  phase  of  the  signal 
source  was  not  steady  enough  over  long  periods.  This  might 
be  another  reason  for  not  having  detected  gravity  wave  winds. 
3esides  these,  from  the  50  m- level  wind  records,  which  did 
not  clearly  show  the  perturbing  winds  of  gravity  waves,  it 
seems  that  a  detailed  theory  of  gravity  waves  propagating  in 
the  turbulent  boundary  layer  is  needed  in  order  to  predict 
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correct  signal  doppler  shifts. 

The  signal  transmitter  has  now  been  put  under  the  con¬ 
stant-phase  control  mode.  This  mode  should  give  a  probing 
signal  with  steadier  phase  over  long  periods.  Three  more 
signal  receivers  are  being  assembled.  The  probe  will  soon 
have  an  array  of  receivers.  These  improvements  of  the  equip' 
ment  should  greatly  increase  the  chance  of  detecting  gravity 
waves. 
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4-5.  Diurnal  variations  of  the  atmospheric  boundary  layer 

The  prolonged  operation  of  the  probe  with  the  receiver 
located  at  the  9.2  Mm  site  has  shown  diurnal  variations  of 
the  atmospheric  boundary  layer.  Various  effective  wind  pro¬ 
files  are  indicated  by  the  diurnal  variation  of  average  signal 
amplitudes.  Usually  no  signal  is  received  from  noon  to  sunset. 
The  receiver  starts  getting  the  signal  around  sunset,  and 
usually  continues  to  have  signal  until  2-4  hours  before  noon 
of  the  next  day.  The  strongest  signals  of  the  day  often  ap¬ 
pear  both  at  the  beginning  and  at  the  end  of  the  receiving 
period.  Fig.  4-5.1  is  a  typical  example  which  shows  diurnal 
variations  of  the  signal  amplitude,  the  wind  speeds  at  3  m- 
and  50  m- levels,  and  the  surface  temperature.  The  example 
includes  two  signal  receiving  periods  on  Aug.  15-17,  1970. 

Four  effective  wind  profiles  at  0700  and  1900  EST  each  day 
were  obtained  by  taking  averages  of  the  regular  upper  air 
soundings  at  Albany,  N,  Y.,  Nantucket  Mass.,  and  Portland, 
Maine.  They  indicate  that  the  night  of  Aug.  16-17  certainly 
had  a  stronger  effective  boundary  shear  than  the  night  before, 
and  that  the  shear  vanishing  height  varies  from  200  m  to  500 
m  on  those  days.  Also  shown  in  Fig.  4-5.1,  are  two  effective 
wind  profiles  from  Boston  pibal  soundings  taken  at  1  am  EST 
each  day.  The  following  is  a  qualitative  discussion  about 
the  typical  effective  wind  profiles  for  different  times  of 
the  day. 

In  the  period  from  noon  to  sunset,  when  no  signal  is 


Fig.  4-5.1  An  example  of  ciurnal  variations 
of  the  atmospheric  boundary  la^er 
on  rtU;.  15-17,  1970 
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received,  the  effective  shear  is  negative  because  of  the 
strong  temperature  lapse  rate. 

At  night,  the  effective  wind  profile  is  possibly  para¬ 
bolic  or  logarithmic  when  the  focusing  factor  is  small,  pre¬ 
sumably  less  than  .5,  and  the  signal  penetrates  up  to  where 
the  shear  almost  vanishes  or  becomes  very  small.  When  the 
focusing  factor  is  large,  presumably  greater  than  .5,  the 
signal  penetrates  only  the  lower  part  of  the  boundary  shear, 
and  the  effective  wind  profile  is  possibly  linear  or  even 
has  a  positive  curvature. 

The  effective  wind  profile,  which  is  responsible  for 
the  strong  signals  at  both  ends  of  receiving  period,  is 
possibly  an  elevated  shear.  As  shown  in  Fig.  4-5.1,  the 
strong  signal  in  the  early  evening  accompanies  the  maximum 
rate  of  decrease  of  the  surface  temperature.  On  the  other 
hand,  the  strong  signal  before  noon  accompanies  the  maximum 
rate  of  increase  of  the  surface  temperature.  In  the  early 
evening,  the  wind  records  show  large  shears  in  the  lowest 
50  m-layer.  But  because  of  the  residual  temperature  lapse 
rate  from  the  daytime  surface  heating,  the  effective  shear 
at  low  levels  is  reduced,  so  that  a  positive  curvature  may 
exist.  Before  noon,  a  large  temperature  lapse  rate  at  low 
levels  is  formed  due  to  the  ground  heating,  and  the  wind  re¬ 
cords  often  indicate  almost  no  shear  in  the  lowest  50  m-layer. 
iis  a  result,  the  effective  profile  before  noon  is  likely 
elevated  shear  type  b  or  c,  which  has  a  zero  or  negative 


surface  shear  (Fig.  3-3.4). 

Diurnal  variations  of  the  signal  amplitude  are  less 
obvious  in  winter  than  summer.  Also  weather  disturbances  can 
obscure  diurnal  variations.  A  summary  of  the  9.2  km  signal 
transmissions  is  shown  in  Fig.  4-5.2,  where  the  focusing 

i 

factors  have  not  been  corrected  for  the  relaxation  damping, 
and  are  smaller  than  the  actual  focusing  factors  by  about  /% 
in  summer  and  30%  in  winter  (section  3-1).  All  the  actual 
focusing  factors  are  smaller  than  unity  except  that  in  the 
early  evening  of  June  12,  1970.  Sunset,  sunrise  and  noon 
appear  to  be  the  natural  dividing  times  between  periods  of 
different  transmission  levels.  The  annual  probability  of 
signal  transmission  to  the  9.2  km  site  is  about  63%  at  night, 
64%  before  noon,  and  27%  in  the  afternoon. 
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Fig.  4-5.2  Signal  amplitudes  receiver  a-  9.2  km  to  the  SE 
of  the  source 
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The  infra-sonic  wave  probe  has  demonstrated  its  capa¬ 
bility  of  measuring  wind  fluctuations  and  turbulence  scales 
in  the  atmospheric  boundary  layer.  The  results  are  generally 
consistent  with  those  of  meteorological  tower  measurements 
and  aircraft  measurements.  The  inferred  wind  fluctuations 
averaged  over  the  entire  ray  path  are  mostly  smaller  than 
the  wind  fluctuations  monitored  at  50  m-level.  "'his  seems 
to  indicate  that  turbulence  is  stronger  at  50  m-levol  than 
at  the  average  signal  propagating  height ,  which  is  estimated 
to  be  about  140  m  to  200  m  from  average  signal  amplitudes 
and  synoptic  meteorological  information.  The  boundary  layer 
eddies  inferred  from  the  probe  data  seem  to  be  more  elongated 
in  horizontal  directions  than  the  eddies  reported  from  tower 
measurements.  The  ratio  of  horizontal  to  vertical  turbulence 
scales  rancres  from  about  1/1  to  If  /I  for  horizontal  turbu¬ 
lence  scales  of  about  200  m  to  1400  m.  The  vertical  turbu¬ 
lence  scale  is  99  +  28  m,  which  is  roughly  57->  of  the  average 
ray  height  and  30-S  of  the  boundary  shear  thickness. 

Diurnal  variations  of  temperature  and  wind  structures 
in  the  atmospheric  boundary  layer  have  been  indicated  by 
average  signal  ampliu'cdes  received  at  a  location  9200  n.  to 
the  Jii  of  the  signal  source.  Jix  fundamental  effective  wind 
profiles  are  used  to  predict  ray  trajectories  and  focusing 
factors  of  the  probing  signal  at  different  times  of  the  day. 

A  shear  vanishing  height  at  about  200  m  to  CC0  n  seems  to  be 
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often  suggested  "by  the  nighttime  probe  data  and  pertinent 
meteorological  information.  Elevated  effective  wind  shears 
are  believed  to  be  responsible  for  strong  signal  amplitudes, 
which  are  frequently  observed  at  both  ends  of  the  daily  signal 
receiving  period,  i.e,,  in  the  early  evening  as  well  as  at 
about  4  to  6  hours  after  sunrise. 

An  array  of  receivers  is  recommended  to  increase  the 
accuracy  of  probe  measurements.  Fy  continuously  measuring 
average  signal  amplitudes  with  several  receivers  located  at 
various  distances  from  the  signal  source,  one  could  investi¬ 
gate  time  variations  of  effective  wind  profiles  in  detail. 

3y  measuring  signal  amplitude  variations  and  phase  shifts  as 
functions  of  the  travel  distance  with  multi-receivers,  one 
cou.ld  test  and  improve  theories  of  signal  propagation  in  a 
random  medium,  fly  studying  spatial  correlations  of  signal 
amplitude  variations  and  phase  shifts,  the  structure  of  atmos¬ 
pheric  turbulence  and  its  effects  on  wave  propagation  could 
be  better  understood. 

The  probe  seems  to  be  ready  to  detect  the  wind  fluc¬ 
tuations  due  to  gravity  waves.  With  an  array  of  receivers 
deployed  at  different  azimuthal  angles  with  respect  to  the 
signal  source,  the  chance  of  receiving  clear  signal  fluctua¬ 
tions  due  to  gravity  wave  winds  would  be  greatly  increased. 

The  probe  could  be  used  in  air  pollution  meteorology. 

The  horizontal  wind  fluctuation  and  vertical  turbulence  scale 
measured  by  the  probe  could  help  predict  the  horizontal  dis- 
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persion  and  vertical  mixing  of  pollutants  in  che  atmospheric 
boundary  layer.  The  shear  vanishing  height,  which  is  indi¬ 
cated  by  average  signal  amplitudes  over  periods  of  hours, 
seems  to  be  associated  with  weakly  turbulent  flows,  and, 
therefore,  could  be  important  to  certain  air  pollution  prob¬ 
lems  (Slade  1969). 
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Appendix  1.  The  probe  data 

This  appendix  is  a  collection  of  tables  listing 
all  available  field  data  and  results  of  analysis.  Examples 
showing  how  to  analyze  field  data  are  given  in  Section  4-1, 
where  the  accuracy  of  these  experimental  results  is  also 


discussed. 
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Table  a— 1.1  Frobe  data  (Kay  £,  1970,  'S  =  4.5  km) 


Period 

big. 

var. 

/7ind 

scale 

_ 

Time 

/V. 

< JF'p 

°$ 

flue. 

hori. 

Vert . 

z 

h 

X 

Xu 

z 

( £3T ) 

(sec) 

(rad) 

(m/s) 

(m) 

(m) 

1705 

153 

.089 

.195 

1706 

160 

.026 

.200 

170  e 

250 

.053 

.425 

Night  188 

aver. 

.056 

.273 

*The  average  signal  amplitude  has  a  focusing  factor  of  1.67. 

A  positive  curvature  profile  with  zero  surface  shear  and  an 
effective  wind  of  7.3  m/s  at  the  300  m- level  (Section  3-3), 
for  example,  could  predict  a  caustic  at.  the  receiver,  however, 
the  simple  ray  theory  can  not  determine  the  exact  value  of 
average  signal  amplitude  near  the  caustic,  'ihe  average  wind 
seems  to  have  blown  along  the  source-receiver  line  in  that 
evening.  The  cross  wind  component  (Vxav)  was,  therefore 
zero.  Then,  according  to  the  model  of  frozen  turbulence 
drifting  with  the  average  wind  (Section  3-5),  the  signal 
fluctuations  were  mainly  caused  by  the  turbulence  near  the 
source  as  well  as  near  the  receiver,  because  of  these 
special  circumstances,  the  field  data  was  not  analyzed  for 
average  horizontal  wind  fluctuations  and  turbulence  scales. 


Period 


3ig.  var 


wind 


Scale 


rj‘xme 


<£/p  <r$ 


flue. 


(ft 


riori , 


Vert, 


L_ 


(HST) 

(sec) 

(rad) 

(m/s) 

(m) 

(m) 

0136 

41 

.074 

.207 

.102 

1058 

101 

0142 

56 

.156 

.276 

.117 

889 

80 

0212 

45 

.074 

.188 

.oee 

1023 

96 

0226 

56 

.156 

.420 

.178 

1070 

99 

0257 

57 

.200 

.344 

.144 

878 

79 

0315 

56 

.156 

.319 

.135 

947 

86 

0333 

46 

.151 

.289 

.134 

901 

84 

0510 

97 

.166 

.353 

.113 

1030 

87 

0514 

96 

.292 

.652 

.210 

1058 

90 

K  ight 
aver. 

61 

.158 

.339 

.135 

982 

89 

107 


‘  able  A- 1.4  Probe  data  (iiug.  16-17,  1970) 


Period 

3ig  . 

var. 

.vine? 

Scale 

Vime 

r 

%/? 

n — 

f  luc . 

Kori . 

Vert. 

°? 

L 

X 

L 

2 

U35  ) 

( sec ) 

(rad) 

(tn/s) 

(m) 

(m) 

2030 

41 

.131 

.126 

.062 

392 

61 

2048 

39 

.082 

.099 

.051 

433 

68 

0047 

32 

.044 

.057 

.031 

436 

70 

0332 

30 

.074 

.062 

.036 

354 

57 

0419 

32 

.094 

.099 

.056 

399 

64 

1042 

26 

.163 

.138 

.104 

206 

66 

1053 

32 

.106 

.163 

.110 

267 

90 

.light 

aver. 

35 

.085 

.089 

.047 

402 

64 

lonv. 

aver. 

29 

.134 

.150 

.107 

236 

78 

able  ii-1.5  Short-period  probe  data  (Nov.  21-22,1970) 


Period 


Time  '?!' 


Sig. 


°P/P 


(EST) 


2155 

2300 

2324 

2351 

0021 

0034 

0232 

0320 

0351 

0430 

0610 

0633 

0740 


( sec) 


69 

50 

50 

45 

56 

50 

38 

56 

45 

50 

45 

41 

45 


.319 

.144 

.138 

.181 

.151 

.099 

.050 

.082 

.082 

.069 

.044 

.050 

.074 


(rad) 


.695 

.289 

.195 

.225 

.262 

.138 

.089 

.138 

.119 

.082 

.074 

.057 

.062 


Wind 

flue. 


0^ 


(m/s) 


.416 

.204 

.136 

.167 

.174 

.097 

.071 

.092 

.089 

.057 

.055 

.044 

.046 


Hori. 

Vert. 

Sc 

Lz 

(m) 

(m) 

0808 

60 

.082 

.106 

.067 

242 

113 

0924 

38 

.057 

.057 

.046 

173 

99 

Night 

aver. 

49 

.114 

.186 

.128 

219 

122 

Conv. 

aver. 

49 

.070 

.082 

.057 

208 

106 

Period 


Sig.  var 


Wind 


Scale 


-ime 

r 

(Tp/p 

0$ 

flue. 

Hori. 

Vert. 

Lx 

LZ 

(  ESi’ ) 

{ sec) 

(rad) 

(m/s) 

(m) 

(m) 

2330 

338 

.41 

1.34 

.41 

1014 

106 

0015 

390 

.28 

.90 

.26 

1170 

105 

0215 

360 

.21 

.58 

.17 

1080 

97 

0230 

255 

.14 

.33 

.12 

765 

91 

0440 

375 

.28 

.76 

.22 

1125 

96 

0530 

300 

.20 

.52 

.17 

900 

95 

0615 

720 

.67 

.97 

.20 

2160 

71 

0740 

900 

.41 

1.52 

.28 

2700 

113 

Night 

aver. 

455 

.325 

.865 

.23 

1364 

97 

Feriod 


Sig.  var 


Wind 


Scale 


Time 

- - 

't 

°P/P 

flue. 

Hori. 

Vert. 

Ov 

B 

B 

(EST) 

( sec) 

(rad) 

(nv's) 

(m) 

(m) 

0916 

41 

.044 

.099 

.080 

200 

152 

0920 

35 

.044 

.082 

.071 

178 

137 

0924 

45 

.069 

.168 

.129 

212 

159 

0928 

38 

.050 

.089 

.074 

181 

136 

0931 

56 

.044 

.138 

.096 

248 

179 

Conv 

aver 

*  43 

.050 

.115 

.090 

204 

153 

'IV' 


N.VI7W 


"M.Ml  WIMi 


1 


111 


Table  a-1.8  Probe  data  (Dec.  14,  1970)* 


Feriod 

3ig. 

var. 

tfind 

Scale 

T  ime 

7: 

0^P 

flue. 

hori. 

Vert. 

<rv 

bx 

L 

2 

( 43T  ) 

(sec) 

(rad) 

(m/s) 

(m) 

(m) 

0413 

ISO 

.370 

1.57 

.605 

360 

233 

0424 

100 

.188 

.482 

.331 

200 

181 

0449 

41 

.238 

.714 

.766 

97 

196 

0510 

75 

.244 

.425 

.337 

150 

150 

Of.  52 

112 

.138 

.539 

.346. 

224 

224 

•light 

aver. 

99 

.260 

.798 

.560 

202 

190 

Conv. 

aver. 

112 

.138 

.539 

.348 

224 

224 

’ Vne  traces  of  signal  records  in  that  night  looked  very  smooth 
(i.e.,  with  little  high  frequency  wiggles),  despite  large 
fluctuations  with  periods  of  a  few  minutes  (Fig.  4-1.3). 

-according  to  wind  records,  there  were  very  little  wind  shears 
and  wind  fluctuations  in  the  lowest  50  m- layer.  Therefore, 
the  signal  variations  could  have  been  caused  by  wind  and 
temperature  fluctuations  associated  with  an  elevated  effec¬ 
tive  wind  shear,  for  which  the  ray  theory  may  not  be  good 
enough  to  predict  amplitude  variations,  -.esides,  no  re¬ 
liable  average  cross  wind  component  (Vxav)  can  be  computed 
from  any  available  weather  information.  rIhe  adopted  model  \ 

with  a  V  of  2.0  m/s  (Table  A-1.10)  is  rather  arbitrary.  \ 

Consequently,  the  analysis  is  only  tentative  and  the  results  j 

may  not  be  meaningful.  £ 


* 

i 
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Table  A— 1.9  Probe  data  (Bee.  19,  1970) 


■I 

Period 

Gig. 

var. 

<^ind. 

Scale 

r 

0^5 

°* 

flue. 

iiori. 

Vert. 

c K 

L 

X 

L 

z 

(EST) 

(sec) 

(rad) 

(m/s) 

(m) 

(m) 

0555 

225 

.363 

.586 

.168 

1153 

80 

0721 

75 

.062 

.326 

.161 

454 

131 

0747 

90 

.062 

.395 

.179 

570 

166 

0025 

56 

.062 

.294 

.168 

316 

136 

0922 

75 

.144 

.457 

.227 

413 

117 

0947 

150 

.464 

.526 

.164 

769 

6S 

1019 

100 

.126 

.502 

.214 

513 

126 

1112 

90 

.126 

.393 

.177 

462 

111 

1147 

90 

.195 

.470 

.213 

462 

97 

1225 

90 

.131 

.432 

.195 

462 

114 

1324 

112 

.126 

.393 

.160 

574 

111 

Night 

aver. 

104 

.147 

.412 

.181 

581 

126 

Conv. 

aver. 

105 

.195 

.453 

.191 

540 

104 

113 


Table  A-1.10  ihe  parabolic  profiles  used  for  the  data 
analysis  in  lables  A-1.2  through  a-1.9 
(V  =*9.2  km) 


r.'ate 

fa 

f 

z 

m 

V 

m 

H 

H 

av 

V 

xav 

(1970) 

(m) 

(m/s) 

(m) 

(m) 

(m/s ) 

(deg) 

6/30 

.40 

.45 

300 

2.2 

201 

141 

9.3 

235 

8/15 

.36 

.39 

290 

2.7 

215 

148 

6.8 

240 

8/16- 

17 

.40 

.43 

372 

3.7 

257 

178 

9.2 

232 

6/17 

.41 

.44 

487 

6.2 

331 

231 

6.4 

260 

11/21- 

22 

.25 

.33 

300 

3.7 

231 

168 

3.0 

270 

11/23 

.29 

.35 

360 

4.5 

275 

194 

3.6 

250 

12/14 

.22 

.30 

375 

5.8 

300 

214 

2_0 

300 

12/19 

.34 

.50 

600 

7.6 

372 

258 

5.1 

234 

Remarks : 

f  =  the  apparent  focusing  factor  inferred  from  the 
a  received  signal  amplitude. 

f  =  the  focusing  factor  estimated  by  correcting  fa  for 
the  dissipations  due  to  relaxation  and  turbulence 
scattering  (Section  3-1). 
zm  =  the  shear  vanishing  height.. 

Vm  =  the  maximum  wind  component. 

H  =  the  penetrating  height. 
riav  =  the  average  ray  height. 

V  =  the  cross  wind  component  at  Ha„. 

XaV  CIV 

ZfL  =  the  prevailing  wind  direction. 
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Table  A-l.ll  Short-period  wind  fluctuations 

(November  21-22,  1970) 


Time 

Period 

Wind  flue. 

at 

7: 

0v 

tfv50 

(EST) 

( sec) 

(m/s) 

(m/s) 

2155 

69 

.42 

1.24 

2300 

50 

.20 

.76 

2324 

50 

.14 

.80 

2351 

45 

.17 

.87 

0021 

56 

.17 

.74 

0034 

50 

.10 

.78 

0232 

38 

.07 

.41 

0320 

56 

.09 

.37 

0351 

45 

.09 

.35 

0430 

50 

.06 

.78 

0610 

45 

.06 

.39 

0633 

41 

.04 

.37 

0740 

45 

.05 

.32 

oeoe 

60 

.07 

.37 

0924 

38 

.05 

.71 

Remarks* 

The  probe  data  are  the  same  as 
listed  in  Table  A-1.5. 

( yv  -  the  rms  wind  fluctuation  inferred 
from  the  probe  data. 

(J^5Q  *  the  rms  wind  fluctuation  monitored 
at  50  m. 


L 


Table  A-l 


Time 


,12  Long-period  wind  fluctuations 
(Nov.  21-22 #  1970) 


Period 


Wind  flue. 


(EST) 


2330 


0015 

0215 

0230 

0440 


0530 


0615 


0740 


(sec) 


(m/s) 


<5*50 

(m/s) 


Remarks : 


( Tv  = 
(Tv50  = 


the  rms  wind  fluctuation  inferred 
from  the  probe  data  of  Table  A-l. 6. 
the  rms  wind  fluctuation  monitored  at 
50  m. 
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Appendix  2.  Ray  characteristics  of  fundamental  effective 
wind  profiles 

THE  NEGATIVE  SHEAR  PROFILE 

For  the  negative  shear  profile, 

V  «  -  Sz,  (1) 

the  ray  angle  is 

Of  =  (.0lo2  +  2  S*/C  )*,  (2) 

which  does  not  vanish  at  any  height.  Therefore  the  ray  will 
never  Lend  back  to  the  Earth's  surface  in  the  negative  shear 
profile. 

THE  PARABOLIC  PROFILE 

For  the  parabolic  profile, 

V  --  S  z  f  1  -  z/(2z  )  ]  =  ( 2V  z/z J  [  1  -  z/(2z  )] , (3) 
0  o  i  m  j  m  mi  ni/ 

the  ray  angle  vanishes  at  the  penetrating  height  H,  which  is 

determined  by 

<*02  *  <2  1  - H  7(2  V]*  °’  (4) 

whence 


where  o[  =(S  z/C)  =(2V/C)  is  the  maximum  value  of 
'm  o  m  o  mo 

C/  for  the  signal  to  bend  back  to  the  Earth's  surface. 

^  o 

The  first  half  trajectory  (from  the  source  to  the 


turning  point)  is 
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and  the  second  half  trajectory  (from  the  turning  point  to  the 
landing  point)  is 


The  entire  ray  trajectory  described  by  Eqs.  (6)  and  (7)  can 
be  expressed  in  a  single  equation,  i.e., 

z  =  zjl  +  <  <Xo/dm)  slnh  <<xmy/z„)  -  cosh  («myAm)]-(8> 

The  travel  distance  is 


■  1"  (U  *0lo/0(m)/a  -0 l0/cxj 

“  tanh  _1  <  •  (9> 
The  average  ray  height  is 


Hav  "  (H [2dz/(  <*o  -  2  W) 

=  2  fl  -  (  a  /c(m)  /  tanh  1  (  OL/oL)  )  .  UO) 
o  'm  o  m  J 

Since  the  derivative 

2^=<2C0/So)/(x-(<Voy2). 

the  focusing  factor  is 

f  =  Y/f$£Kf 

=  {orm[  1  -  (Otar's  A2ao)}h{in  (U *oj/ctm)/U-a[/ctm))f 


which  is  a  function  of  either  (  oL/oO 

o 

or  OiY  /(2z  )  =  (y/(2C  )*  1  V J5  /  z. 
m  m  v  o  J  m  m 

With  given  Y  and  CQ,  the  f  -  contours  on  the  log  Vm  vs.  log 

zm  diagram  are  straight  lines  (Fig.  3-3.2). 

The  ray  heights  as  fractions  of  zm  are  functions  of 

(  oL/of  )  and,  hence,  are  also  functions  of  the  focusing 
o  *m 

factor  f  (Fig.  3-3.3). 

THE  LOGARITHMIC  PROFILE 

For  the  logarithmic  profile 
Ve=Vt  ln(Z/z0), 


(12) 


i 


the  ray  angle  can  be  normalized 


=  (g  2  -  In  (z/ze)  ) 


(13) 


where  G  -  oi/olj,  , 

Go  =  V°(e  ' 

and  0(g  =  (2  Vj/C^  . 

The  penetrating  height  is 
H  =  exp  (Gq2). 

The  trajectory  from  the  source  to  the  turning  point  is 


(14) 


=  U/Cfc)  f 

A z_ 


(dz/G) 


=  (-2zQ/0jg)  exp  (Go2)  /  exp  (-  G2)dG 


=  (  (7J5  ZQ/ Ojj)  exp  (GQ2)](erf  (GQ)  **  erf  (G)  )  ,  (15) 

where  erf  (G  )  =  ( 2/i h  fQ  °  exp  (-  u2)  du. 

*o 

The  trajectory  from  the  turning  point  to  the  landing  point 
is  similarly  obtained  as 

3g  2 

y  —  f  (TT  zo/0(j)  exp  (G  )]  (erf  (G  )  +  erf  (G) ) .  (16) 


The  travel  distance  is 

2 

Y  =  (2  7T  z  / CL )  exp  (G  )  erf  (G  ). 

o  •*  o  o 


(17) 


The  average  ray  height  is 
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H  =  z  exp  (G  2)  erf  (2^  G  )  /  (  2^  erf  (G  )]  (18) 

av  o  o  o'  o  * 

Since  the  derivative 

=  <“  Go  “P  O  erf  <Go>)  • 

the  focusing  factor  is 


(19) 


In  a  typical  case,  2  m/s,  0^^331  m/s.  Then  the 
normalizing  angle 
Oljl  =  .11  rad. 

With  an  initial  ray  angle  0/o*S?  -15  rad»  which  is  normalized 
as  G^=»  1.4,  the  error  function 

O 

erf  (GQ)  1 

(the  error  being  less  than  5%,  see,  e.g.,  Abramowitz  et  al 


1964,  p.  311).  As  a  result,  the  following  approximate  for¬ 
mulas  are  applicable  for  Y  j=j  250  m. 


Y  =  2  IT  H/  Oil , 

H  =  H/2^, 

av  ^ 

f  ~  1/^2  In  (  0(£  Y/  ( 2  77^  ZQ))J 


(20) 

(21) 

(22) 


THE  LINEAR  PROFILE 

For  the  linear  profile 
Ve  =  sz. 


(23) 
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the  penetrating  height  is 

H  =  Cc  tfQ2/(2S).  (24) 

The  trajectory  is 

y  =  (2  Cq/S)  IH^  +  (H  -  z)  J 

or  (y  -  C  0 LAO2  =  (2  C  /S)  (  C  o/ 2/(2S)  -  z  ]  .  (25) 

The  travel  distance  is 

Y  =  2  C  OL/S.  (26) 

o  o 

The  average  ray  height  is 

Hav  “  co^o2/(3s)  (27) 

=  ( 2/3  )H 

The  focusing  factor  is 

f  =  1,  (28) 

for  which  the  second  order  term  can  he  determined  by  first 
integrating  the  ray  trajectory  without  shallow  angle  approxi 
mation  and  then  expanding  the  result  in  small  ray  angles. 
Thus  one  gets 

f  =  1  -  5S2  Y2  /  (32  Co2),  (29) 

which  decreases  very  slowly  as  Y  increases. 

THE  POSITIVE  CURVATURE  PROFILES 

The  ray  characteristics  of  the  positive  curvature 
profile  type  a 

Ve  =  Sc  z  (  1  ♦  z/(2zp)  )  ,  (30) 
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can  be  obtained  from  the  corresponding  formulas  of  the  para- 

Jj 

bolic  profile  by  putting  z„  =' -  z  and  N  =  (-1)  N  , 

m  p  ^  p 

where  o/p  =  (SD  zp /ZQ)h. 

For  the  positive  curvature  profile  type  b 

Ve  =  SQ  z2/(2zp),  (31) 

one  gets  the  penetrating  height 


H  =  Qk  w 

where  ofp  =  (so  zp/co)  ; 
the  ray  trajectory 

y  =  (zp/o(p)  (“W/ 2  +  cos-1  (z/H))  , 
or  z  =  H  sin  (  0(pY/zp)  ? 

the  travel  distance 

Y  =7Tz  /o'  * 

P  P 

the  average  ray  height 

Hav  =  (2/1r)  H 
the  focusing  factor 
f  =  OO  # 


(32) 


(33) 


(35) 


(36) 


which  means  that  a  caustic  is  predicted  at  7T  z  /  qfD  from 
the  source. 

For  the  positive  curvature  profile  type  c 


ve  =  V  f  -  1  +  z  /<2zp> )  , 

the  penetrating  height  is 


H  =  Z 


,{  i  +(i  ♦  (or0/ <xp)2)H]  . 


(37) 


(38) 
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The  first  half  trajectory  is 


ana  the  second  half  is 


The  entire  trajectory  can  be  expressed  as 

z  =  z  [  1  +  (  Of/  0(  )  sin  (  Of  y/z  )  -  cos  (Of  y/z  )l(4l) 
pi  op  PP  p  p  J, 

The  travel  distance  is 

Y  =  (2  zp/tfp)(TT/2  ♦  cot”1  (P^/0/p)).  (42) 

The  average  ray  height  is 

Hav  *  zp{  1  +  <  ao/%)/  (lr/2  *  cot'1  (Q,</«p>)j  •  (43) 

The  focusing  factor  is 

£  =  «y*o)!5 11  +  moM)2]'‘  t7172  * cot'1  (‘V<vf'(44) 

which  predicts  two  caustics  at  the  distances  of  TT  z  /tf  and 

P  P 

2  TT  z  /cl  from  the  source.  All  the  rays  land  between  1 T  z  /#  , 
P  P  P  P 

&2  TT  Zp/Qp#  The  theory  also  predicts  a  skip  distance  of 

TTVV  since  no  ray  is  received  within  this  distance. 
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THE  ELEVATED  SHEAR  PROFILES 

The  ray  characteristics  for  elevated  shears  can  be 
evaluated  numerically  (Wesson  1970) ,  and  the  results  are 
schematically  shown  in  Fig.  3-3.4. 

The  interesting  feature  of  elevated  shear  types  b  and 
c  is  the  skip  distance  Y^  where  a  caustic  is  also  pre¬ 
dicted.  can  be  estimated  by  using  a  two-layer  model. 

The  lower  layer  of  a  thickness  of  h^  has  a  negative  shear 
-S^  and  the  upper  layer  has  a  positive  shear  S2.  Then  the 
horinontal  travel  distance  in  the  lower  layer  is 

(2  C</S1)  (  (  0to2  *  2Sl  hj/co)  -  a'J  . 
and  that  in  the  upper  layer  is  from  Bq.  (26) 

(2  C/S.,)  (  O/J  *  2SX  VCo,!S- 
The  total  travel  distance  is  the  sum  of  these,  i.e.. 


Y  = 


2C, 


u.%  (oi2  -  ajaJ 


V  -o 


.  (45) 


3y  putting  the  derivative  QL.  to  vanish  and  solvincr  for  the 

30lo 

value  of  C (Q,  one  gets 

‘°loW  *  (2S1  VCo)l5  '  ((1  +  VS2>2  ~l)H  <46) 


Substitution  of  £q.  (46)  into  Eq.  (45)  gives 


1  2  1 

Ymini  =  2  ^W^l*  ((l+Sj/S^  ”1  X  •  (We  c>  <47) 


r7«'  ty’tt’"' 
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l 


which  is  the  skip  distance  for  elevated  shear  type  c.  The 
skip  distance  for  elevated  shear  type  b  is  obtained  from  Eq. 
(47)  by  setting  S1  to  zero,  i.e., 

,5S 


Ymini  =  4  <co  VS2J 


(type  b) 


(48) 


? 

1 


& 
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Appendix  3.  Parabolic  profile  fluctuations 
According  to  the  linear  theory,  the  time  rate  of  frac¬ 
tional  change  of  amplitude  is 

fa  Un  p> 

=  s°  It  (ln  So)  4:  tln  p)  +  zm  ft  <ln  zm>  ~<1»  P>.  U) 

And  che  fractional  '-'oppler  is 

£r>  =  A  u)/u) 


=  <ln  so)  I-so  +  ~  <ln  2m>  Ezm' 
where,  according  to  Eq.  (4)  in  Section  3-1, 


dSo  =  (2s0/c02)  f  dz 


(2a) 


ane  :  zm  =  <2zn/co2'  fH  dz  ^i§/0()*  (2b) 

c  ®  m 

Here  the  effect  of  penetrating  height  variations  on  the 

doppler  is  negligible  to  the  first  order. 

To  calculate  —  (In  p)  (M  stands  for  Sr  or  zm),  one 
dM  '  1,1 

should  consider  that  both  source  and  receiver  are  fixed  in 
position,  and,  therefore,  the  derivative  includes  wo  .ertns, 
i.e. , 

1.  =  _J_  (3) 

a.M  »M  5 0to  ' 

The  first  term  is  due  c  .he  explicit  dependence  of  p  .-n 
The  second  term  arises  'ecause  of  he  adjustment  of  :  he 
initial  rav  an  le,  0fo,  to  insure  a  constant  horizontal  travel 
distance,  ,  while  M  is  varvin-.  A  constant  *  reans  that 


-w 1 1  HI »  m  lyPLiPniiiw^^iwsiUP  pwf 


3  «*».*.*  ^T-^fPSIy 
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Ik  «  *  M  *<*°  -  °* 

whence 

3c/o  _  y_-  ,  _a  y 

*dM  "  IK  7  3c 10  * 

Now  from  Eq.  (7)  in  Section  3-2,  one  obtains 


(4) 


(5) 


r?  v  _  -1 

dK  ( 1°  P )  **  o 


3  v 
ddG 


=  zif/jd.v_ 

2  lAaii  ao(0 


3o(c  3*Y 
3M  3(*02 


Mg 

3m' 


0foJ  (6) 


After  reing  the  algebra  of  Eqs.  (6)  anri  (2),  one  .;ets 

-so  ?f~  <ln  P}  =  zm  (ln  p) 

=  (^)[(S0Y/(2  Co0fo))(l  +  (0(/Olm}2j-  l](7a) 
~So  co)J^(so-/ ( 2  ccAo)}(l  +  (°(yCtfm)  }-lj  (7b) 

Dzm  =(3  zm  0^/(2  C0)j[(s0V(2  Cc<yo;J  (l  -  (0^/0^)  2/3j-l](  7c) 


Substitution  of  Eq.  (7)  into  Eqs.  (1)  anr*  (2)  gives  Eqs.  (1) 
ana  (2)  in  Section  3-4. 
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one  can  chance  the  integration  variables  of  Kq.  (4)  into 
<o  =  +  ^2)/2 

=  the  center-cf-mass  coordinate, 


=  the  separation  coordinate, 

and 

(&<p)2  =(0Q2  (4\ve)2/C4J f  dy0  f  dy'. 

o  -Y 

.  exp(-  (x‘/ax)2  -  (y'/ay)2  -  (z'/a2)2  ]  (6) 

The  integration  limits  for  y  can  be  set  to  +  oo  ,  since  Y  i.s 
much  greater  than  ay.  To  the  first  order,  the  integral  can 
be  evaluated  at  x1  =  z1  =0,  assuming  that  the  deviations  of 
the  actual  ray  trajectories  from  the  average  trajectory  are 
small  in  comparison  with  ax  and  a2.  Then  Eq.  (6)  can  be 
evaluated  as 

(A  <p)2  =  (77^  CO2  ay  Y/34)  (A  ve)2.  (7) 

With  the  y  -  coordinate  (Fig.  3-1.1)  as  the  integration 
variable,  the  signal  amplitude  is  given  by 

’■(’’"•'H&Saf. 

where  Of©  the  initial  vertical  spreading  angle, 

1f0  =  If  -  p0 

=  the  initial  spreading  angle  in  the  x-direction. 
Then  the  fractional  amplitude  fluctuation  is 


(9) 
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The  average  value  of  the  derivative  M-  can  be  ob,aine,  from 

the  trajectories,  which  are  evaluated  for  fundamental  effective 

wind  profiles  In  section  3-3.  -she  averaoe  M-  is  equal  Y 

,  ,  Wo 

the  same  as  m  a  uniform  medium.  The  randomly  fluctuating 

derivatives,  A  and  1  „ 

\dcioJ  *{iroj'  are  now  to  136  derived.  For 

'  °ne  needs  Z-as  an  integral  in  y,  which  is  obtained 

from  (2)  in  Section  3-1  and  the  shallow  angle  approximation 


2  -/  ^  ((c  sino/  +  V2)/<c  cosoU  V  )] 

yo  y  J 

=  f  *y (<*-  (€  +  vy  +  Vy  +  vz)/c J  , 


)/c)  , 


(10) 


whence 


3<tfo  'I  Hmo  -  dk  16  *  vy  +  vy  +  vt>/c  j 


(id 


The  elevation  angle  d  Is  from  Eqs.  (2)  and  (3)  in  section  3-1 


«=Ct°  +  £Y 

mcte  *fj  flVk~  ^  ~  *’  £)//<:  *  Vj 

=  do-(  1/5)  fV  dy/’H  +1L;  +  ?i) 

/  \  3Z  9  Z  32  /  ' 


(12) 


whence 


~  =  1  “  (1/C) 


/*'  dJJL  +  ?Vy 

/  9  z  +  3z, 


(13) 
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.(5vV<4  C2  Y2) J  /"  dyc  /*  cv*  (y^2  -  y'2/4) 

b  - Y 

*  (v2  -  2y0  Y  +  yQ2  -  y,2/4)  • 

jj£n  +  ^72) exP  f  -  <x'/ax)2  -(y'/av)2  -  (z'/a. 


)2  -  (z>/az)2J  (18) 


Here  the  center-of-mass  coordinate,  rQ,  the  separation  coor¬ 
dinate,  r',  and  the  oaussian  correlation  function  have  been 
used.  By  noting  that  Y  »  ay,  and  assuming  small  deviations 
of  actual  ray  trajectories  from  the  average  trajectory,  one 
can  evaluate  the  integrations  of  Eq.  (18)  and  obtain 


(*P)2/P2  =((Ave)2  a y/c2)- 

•  [  (7T2  \,3/10)/at  +  (V*X3/15)  /  ( az2  ax2) 

+  (TT^'VlO)  /  ax4  J 


(19) 


Similarly,  one  can  derive  the  mean  square  fractional 
amplitude  fluctuation  for  t-he  parabolic  profile  as 

(Ap)2/ p2  =  |  (4ve)2  ay/^2J* 


*^.<1  /  az4  +  K2/(az2  ax2)  +  K3  /  ax^  J  , 


(20) 


where 


%  =  ^3  7 T5  V  (z^/fc^)2/  (4  sinh  2(Y*)J  J  • 

•  |sinh  (2  Y*)  /  (2  Y*)  -  1  -  (2/3)  (V*)2  exp  (ay*/2 )2J 
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IC2  =  ( 4  TT*5  exp  (ay*/4)2  (zm/tfm)3  /  sinh(Y*)J. 

•(2  +  cosh  (Y* )  -  (3  /  y*)  sinh  (Y*)J  , 

K3  «  IT*  Y3  /  10, 

y*  =  *  <4  /  zn' 

ay  =  ay  P^m  ^  zw 

<4=  (2Vm/0o))S=  (30zm/C0)!S, 
exp  (ay*  /  2)2~  exp  (ay*  /  4)2~  1 


REFERENCES 


Abramowitz,  Milton  and  Irene  A.  Stegun  1964 
Handbook  of  Mathematical  Tables 

Use?iesa55°(?964)reau  "  Stanaards  APPlied  Mathematics 

Batchelor,  G.  K.  1957 

Wave  scattering  due  to  turbulence 

Spp^°409~430  NaVal  Kydr°dynamics  (F*  s*  Gherman,  Ed.,  1957) 
Bunker,  Andrew  F.  1955 

T“f?ulf?ce  and  scaring  stresses  measured  over  the  North 
tA  iC  Ccean  hY  an  airplane-acceleration  technioue 
J.  Meteor.  Vol.  12  (Oct.  1955)  pp.  445-455  tecnnique 

Byzova,  N  L.,  v.  N.  Ivanov  and  S.  A.  Morozov  1965 

characteristics  of  the  wind  velocity  and  temperature 
PrioCt^?ti0n%in  the  a tnio spheric  boundary  layer 

^turbulence  an /r* Jnternational  colloquim  on  atmospheric 
?nd  radl°  wave  propagation  (Moscow,  1965) 
publishing  house  "NAUKA"  Moscow  (1967)  pp.  76-92 

Chernov,  Lev  A.  1960 

Wave  propagation  in  a  random  medium 
Dover  Pub.,  Inc.  New  York  (1960) 

Cole,  John  E.  Ill  and  Richard  A.  Dobbins  1970 
Propagation  of  sound  through  atmospheric  foa 
J.  the  Atmos.  Sci.  Vol.  27  (1970)  pp.  426-434 

Daniels,  Fred  B.  1959 

70if^redUSin9  line  ">icroPh°ne  for  frequencies  below  1  cps 
J.  Acous.  soc.  Am.  Vol.  31,  Ho.  4  (April  1959)  pp.  529-531 

Gardner,  Floyd  Martin  1967 
Phase lock  techniques 
New  York  Wiley  (1967) 

Haurwitz,  Bernhard  1941 
Dynamic  Meteorology 
McGraw-Hill  Book  Company,  Inc.  (1941) 

Hayes,  Wallace  D.  1970 
.Kinematic  wave  theory 

Proc.  Roy,  Soc.  Lond.  Vol.  A320  (1970)  pp.  209-226 

Huschke,  Ralph  2.  1959 
Glossary  of  Meteorology 
Am.  Meteor.  Soc,,  Boston,  Mass.  (1959) 


135 


Knesisr,  H.  C.  1965 

Relaxation  processes  in  gases 

“Physical  Acoustics'*  (ri.  P.  Mason,  Ed.),  Academic  Press 
Inc.,  New  York  Vo 1.  2A  (1965)  pp.  133-202 

Lettau,  Heinz  H.  and  Ben  Davidson  1957 
Exploring  the  atmosphere's  first  mile 
Pergamon  Press  (1957)  Two  volumes 

(Aircraft  measurements  are  described  on  pp.  267-275, 
471,  and  495-496) 

Lighthill,  M.  J.  1965 
Group  velocity 

J.  Inst.  Math.  Appl.  Vol.  1,  (1965)  pp.  1-28 

Lumley,  John  L.  and  Hans  A.  Panofsky  1964 
The  structure  of  atmospheric  turbulence 
Interscience  Publishers  (1964) 

Madden,  T.  R.  and  J.  F.  Claerbout  1968 

Jet  stream  associated  gravity  waves  and  implications 
concerning  jet  stream  stability 
Symposium  proceedings  on  "acoustic-gravity  waves  in  the 
atmosphere"  at  Boulder,  Colorado  on  15-17  July  1968 
(edt,  T.  M.  Georges)  pp.  121-134 

Morse,  Philip  M.  194S 
Vibration  and  sound 

McGraw-Hill  Book  Company,  Inc.,  2nd  edition  (1948) 

Morse,  Philip  M.  anc  K.  Uno  Ingard  1968 
Theoretical  acoustics 
McGraw-Hill  Book  Co.  (1968) 

Cbuchow,  A.  M.  1953 

tiber  den  Einfluss  schwacher  Inhomogenitcften  der  Atmos¬ 
phere  auf  die  Schall-  und  Lichtausbreitung 
"Sammelband  zur  statistischen  Theorie  der  Turbulenz" 
(edited  and  translated  by  Herbert  Goering)  Akademie- 
Verlag  Berlin  (1958)  pp.  157-171 

Olson,  Harry  F.  1947 

Elements  of  acoustic  engineering 

D.  Van  Nostrand  Co.,  Inc.  2nd  Ed.  (1947) 

Panofsky,  H.  A.  and  I.  A.  Singer  1965 
Vertical  structure  of  turbulence 

Quart.  J.  Roy,  Meteor.  Soc.  Vol.  91  (1965)  pp.  339-344 

Pielke,  R. A.  and  H.  A.  Panofsky  1970 

Turbulence  characteristics  along  several  towers 
Boundary-Layer  Meteor.  Vol.  1  (1970)  pp.  115-130 


136 


Piercy,  J.  E.  1969 

Role  of  the  vibrational  relaxation  of  nitrogen  in  the 
absorption  of  sound  in  air 
J.  Acous.  Soc.  Am.  Vol.  46  (1969)  pp.  602-604 

Slade,  David  H.  1969 

Low  turbulence  flow  in  the  planetary  boundary  layer  and 
its  relation  to  certain  air  pollution  problems 
J.  Appl.  Meteor.,  £  (Aug.  1969)  pp.  514-522 

Tatarski,  V.  I.  1961 

Wave  propagation  in  a  turbulent  medium 
Dover  Pub.,  Inc.  New  York  (1961) 

Taylor,  G.  I.  1935 

Statistical  theory  of  turbulence 

Proc.  Roy.  Soc.  Lond.  Vol.  A151  (1935)  pp.  421-476 

Thuillier,  R.  H.  and  U.  O.  Lappe  1964 

Wind  and  temperature  profile  characteristics  from 
observations  on  a  1400  ft  tower 
J.  Appl.  Meteor.  Vol.  3  (June  1964)  pp.  299-306 

Wesson,  Robert  L.  1970 

A  time  Integration  method  for  computation  of  the  intensi¬ 
ties  of  seismic  waves 

Bull.  3eis.  Soc.  Am.  Vol.  60,  No. 2  (1970)  pp.  307-316 


ACKh  CWLEDG  EMSNTS 


I  am  most  indebted  to  Frofessor  T.  r.  Madden  for  his 
guidance  and  encouragement  throughout  the  course  of  this 
thesis,  he  designed  all  of  the  electronic  circuits,  and 
shared  many  of  the  exploratory  field  trips  with  me.  i  would 
also  like  to  acknowledge  very  helpful  discussions  with  K. 

Aki'  JOhn  3r0wn'  L*  Dean,  Tom  De  Fazio,  Dave  Fitterman, 
D.  F.  Keilv,  E.  L.  Mollo-Christensen,  A.  D.  Fierce,  F. 

Sanders,  and  Ron  Ward. 

i-ly  sincere  thanks  also  go  to  3am  Hendryx,  George  r-umont, 
and  many  others  who  have  kindly  assisted  in  carrying  out  the 
experiments. 

I  express  my  affectionate  gratituoe  for  my  wife’s 
understanding ,  sacrifice,  and  skillful  typing,  and  to  my 
late  father-in-law,  Mr.  Yen-chin?  Lin,  without  whose  Genero¬ 
sity  my  study  at  M.I.l.  would  never  have  been  possible. 

General  support  was  received  from  the  U.s.  Army  Re¬ 
search  Office  Contract  s o,  DA-31-124-ARC-431. 


Trrr- 


rwff'wjlP 


138 


BIOGRAPHICAL  NcTE 

The  author  received  his  undergraduate  education  at 
National  Taiwan  University,  China,  and  obtained  a  K.S.  in 
Physics  in  1961.  He  attended  graduate  school  at  National 
Tsing-Hua  University,  Taiwan,  China  from  1961  to  1963  (to.s.) 
as  well  as  at  MIT,  U.S.A.  from  1965  to  1972  (Ph.P.). 

The  author  was  married  to  Shu-Mei  Lin  in  1964.  They 
have  a  daughter  named  Iris  Ting-Lan  Chung. 


